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rovides  a  very  smooch  and  clean  container  for  the  levitated  drop,  and  because 
the  shape  oscillation  is  excited  acoustically,  without  contact  with  a  solid 
surface,  the  drop  may  survive  in  the  superheated  state  for  a  long  enough  time 
for  the  experiments  to  be  conducted. 

A  simplified  model  was  developed  for  describing  the  free  oscillation  and 
forced  oscillation  of  a  drop  freely  suspended  in  a  host.  The  results  of  the 
simplified  model  agree  well  with  that  of  previous  theories. 

The  accuracy  of  this  technique  was  tested  by  measuring  the  interfacial 
tension  between  water  and  each  of  two  common  liquids.  Experiments  were  con¬ 
ducted  for  measuring  the  interfacial  tension  between  water  and  each  of  the 
three  superheated  liquids;  butane,  isobutane,  and  propane,  at  room  temperature 
and  atmospheric  pressure.  Data  for  the  compressibility  of  these  superheated 
liquids  were  also  obtained  using  an  established  acoustic  levitation  technique. 
Some  other  possible  applications  of  the  drop  deformation  technique  are 
mentioned. 
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Chapter  1 


INTRODUCTION 

In  various  fields  like  physical  chemistry,  chemical  engineering,  and 
fluid  mechanics,  surface  and  interfacial  tension  have  long  been 
interesting  subjects.  They  have  proven  to  be  of  fundamental  importance 
in  studies  of  the  properties  of  solutions,  intermolecular  forces, 
molecular  orientation  on  liquid  surfaces,  colloid  phenomena, 
capillarity,  phase  transitions,  and  numerous  aspects  of  theoretical  and 
industrial  chemistry.  Various  methods  of  measurement  have  been 
developed  and  numerous  results  of  the  surface  and  interfacial  tensions 
of  different  liquids  have  been  reported  in  the  literature.  On  the  other 
hand,  for  the  situations  where  one  liquid  is  in  metastable  state,  say 
superheated,  the  traditional  methods  for  measuring  interfacial  tension 
are  unsuitable  and  data  are  unavailable.  Also,  none  of  the  existing 
theories  has  been  proved  universally  applicable  in  predicting 
interfacial  tensions  in  various  circumstances.  In  this  work  we  study  an 
experimental  technique  for  measuring  interfacial  tension  and  apply  it  to 
determine  the  interfacial  tensions  between  water  and  each  of  three 
superheated  liquids.  Another  bulk  property,  the  compressibility,  of  the 
three  superheated  liquids  is  also  studied  using  a  similar  technique. 


1.1  Interfacial  Tension 


When  a  liquid  touches  another  fluid,  a  boundary  forms  and  a  transitional 
region  is  established.  Generally,  we  call  the  contact  region  an 
interface  if  the  other  fluid  is  a  liquid,  or  a  surface  if  the  other 
component  is  a  gas.  Interfacial  and  surface  tension  are  natural 
consequences  of  the  differential  attraction  force  across  the 
transitional  region.  For  example,  let  us  first  consider  the  simpler 
one,  surface  tension.  In  each  bulk  phase,  the  molecules  are,  on  the 
average,  subject  to  a  spherically  symmetric  attractive  force  field.  The 
force  is  much  stronger  in  the  liquid  phase  than  it  is  in  the  vapor  or 
gas  phase.  Because  of  the  difference  in  forces,  the  molecules  in  the 
interfacial  region  are  subject  to  an  asymmetric  intermolecular 
attraction.  Therefore,  work  must  be  done  against  the  attractive  forces 
within  the  interior  of  the  liquid  to  bring  molecules  to  the  surface 
region.  This  work  is  called  the  free  surface  energy  and  can  be  defined 
as  the  work  required  to  increase  a  unit  surface  area  at  constant 
temperature.  The  commonly  used  dimension  of  free  surface  energy  or 
surface  tension  is  ergs  per  square  centimeter  or,  in  SI  units,  newtons 
per  meter.  The  potential  energy  decreases  and  the  stability  increases 
as  the  surface  area  decreases.  This  implies  the  tendency  for  the 
surface  to  attain  a  minimum  area,-  this  is  essentially  where  the  term 
surface  "tension"  came  from.  The  transitional  region  is  estimated  to  be 
about  two  molecules  thick  [1].  Compared  with  most  interesting 
dimensions  in  fluid  mechanics  and  certainly  in  this  v/ork,  this  thickness 
is  so  small  that  the  transitional  region  can  be  treated  as  a 
discontinuity. 


The  origin  of  interfacial  tension  is  the  same  as  that  of  surface 
tension,  except  that  it  is  concerned  with  both  cohesional  forces  between 
similar  molecular  species,  and  adhesional  forces  between  two  different 
species.  The  magnitude  of  the  cohesion  and  adhesion  could  be 
comparable.  It  is  reasonable  to  expect  that  the  interfacial  tension  o 

AB 

between  liquids  A  and  B  will  be  no  more  than  0  +  0  ,  where  0  and  0 

A  B  A  B 

are  the  surface  tensions  against  their  vapors  of  liquids  A  and  B.  0. 
will  equal  0.  +  o_  if  adhesion  is  negligible.  Therefore,  in  general, 
the  increase  of  similarity  in  structure  between  two  liquid  species  i 
reflected  in  a  smaller  value  of  interfacial  tension  and  a  greater 
tendency  for  mutual  molecular  dispersion,  that  is,  miscibility.  The 
interfacial  tension  must  be  positive  for  the  stable  existence  of  an 
interface. 

1.2  Traditional  Methods  for  Interfacial  Tension  Measurements 

The  methods  of  measuring  surface  and  interfacial  tension  may  be 
generally  divided  into  two  classes,  dynamic  and  static.  In  the  dynamic 
methods,  the  interface  periodically  expands  and  contracts,  and  the 
interfacial  tension  operates  to  restore  the  interface  to  a  minimum  area. 
By  measuring  the  parameters  of  the  periodic  motion,  such  as  the 
frequency,  the  interfacial  tension  can  be  calculated.  The  ones 
mentioned  in  the  literature  involve  the  method  of  oscillating  jets, 
oscillating  drops,  and  ripples  [1].  In  the  static  methods,  quiescent 
interfaces  are  involved,  and  the  interfacial  tension  plays  the  role  of 
balancing  other  forces  on  the  interface  due  to  either  the  gravitational 


field  or  the  pressure  differential  across  the  interface.  The 
capillary-height  method,  the  ring  detachment  method,  and  the  drop-weight 
method  are  three  examples  [1] . 

In  general,  dynamic  methods  do  not  give  results  as  reproducible  and 
consistent  as  static  methods.  Also,  the  results  obtained  by  dynamic 
methods  do  not  generally  agree  with  those  obtained  by  static  methods. 
This  discrepancy  is  especially  significant  when  either  of  the  two 
liquids  is  a  solution  of  surface  active  solutes.  When  a  solution  is  in 
equilibrium,  the  composition  in  the  interfacial  region  is  different  from 
that  in  the  bulk  phase.  The  true  interfacial  tension  can  be  obtained 
only  when  all  involved  liquids  are  in  equilibrium.  Therefore,  if  the 
time  required  by  the  diffusion  of  the  solute  molecules  ‘•o  establish  the 
state  of  equilibrium  is  significant  compared  to  the  time  scale  of 
measurement,  dynamic  methods  can  not  give  the  true  interfac.ial  tension 
[2]  [1]. 

1.3  Superheated  Liquids 

At  atmospheric  pressure,  100°C  is  referred  to  as  the  normal  boiling 
point  of  water.  It  is  also  known  that,  under  the  same  pressure,  there 
are  circumstances  in  which  water  can  be  heated  beyond  100°C  without 
boiling;  we  call  water  in  this  state  "superheated".  A  superheated 
liquid  is  one  particular  instance  of  a  metastable  state  of  matter.  In 
thermodynamics,  the  stable  states  of  matter  are  separated  by  equilibrium 
lines  in  pressure-temperature  coordinates.  Matter  that  starts  from  a 
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stable  state  and  transgresses  an  equilibrium  line  without  phase 
conversion  enters  a  metastable  state.  How  far  the  state  deviates  from 
equilibrium  is  characterized  by  the  depth  of  penetration--the 
corresponding  temperature  or  pressure  difference.  A  metastable  state 
can  survive  for  a  very  long  time  even  under  large  variations  in  the 
thermodynamic  parameters.  On  the  other  hand,  the  metastable  state  may 
revert  to  a  stable  phase  if  that  phase  (like  vapor  bubbles  in  the  above 
example  of  superheated  water)  is  introduced  in  some  way  into  the 
metastable  state. 

Generally,  phase  transitions  can  ha  either  of  two  types.  The  first 
type  is  homogeneous  nucleation.  This  happens  if  the  nucleus  of  the  new 
phase  arises  spontaneously  due  to  thermal  fluctuations  and 
intermolecular  interactions.  But  usually  the  phase  transition  starts  at 
the  walls  of  the  container  or  on  foreign  inclusions.  This  type  of 
process,  called  heterogeneous  nucleation,  requires  a  much  lower  degree 
of  superheat  [3].  Therefore,  one  expects  that  a  difficulty  in 
experimental  studies  of  superheated  liquids  is  the  prevention  of 
heterogeneous  nucleation.  Metastable  liquids  are  of  interest  in  fields 
like  molecular  physics  and  in  problems  such  as  the  study  of  the  initial 
stage  of  phase  transition,  where  the  interfacial  or  surface  tension  is  a 
very  important  parameter.  The  vaporization  of  superheated  liquids  may 
be  explosive  in  character.  One  of  the  motivations  for  this  work  is  that 
the  unintentional  contact  of  a  cold  liquid,  like  liquified  natural  gas 
(LNG) ,  with  a  much  warmer  stable  liquid,  like  water,  may  cause  dramatic 
vapor  explosions  [4] .  The  interfacial  tension  may  be  needed  in  modeling 


this  phenomenon  if  the  superheat  mechanism  is  assumed  to  be  the 
initiator  of  the  physical  explosion. 

1.4  Acoustic  Levitation 

To  measure  the  interfacial  tension  between  water  and  a  superheated 
liquid,  most  of  the  commonly  used  techniques  have  to  be  ruled  out 
because  of  the  involvement  of  solid  surfaces  which  would  trigger  the 
vaporization  of  the  superheated  liquid.  The  idea  of  oscillating  drops 
seems  to  be  the  natural  choice.  If  we  can  put  a  drop  of  a  superheated 
liquid  (say  a  hydrocarbon)  into  an  immiscible  "host"  liquid  (say  water), 
then  the  host  will  act  as  a  very  smooth  and  clean  container  for  the 
drop,  thereby  reducing  the  probability  of  the  presence  of  solid 
heterogeneous  nucleation  sites.  If  a  hydrocarbon  drop  oscillates  in 
water,  then  the  interfacial  tension  which  provides  the  restoring  force 
in  the  interfacial  disturbances  could  be  inferred  from  some  measurements 
of  the  oscillation. 

The  method  which  we  use  to  introduce  a  superheated  hydrocarbon  drop 
into  water  will  be  described  in  chapter  4.  But,  how  do  we  excite  the 
drop  into  oscillation  and  take  measurement  of  the  oscillation?  Acoustic 
levitation  is  a  very  appropriate  technique  for  this  situation  [5].  A 
strong  acoustic  standing  wave  can  be  established  in  a  host  liquid,  water 
in  this  work,  and  produce  an  acoustic  force  on  an  immersed  drop.  By 
adjusting  this  force  to  balance  the  gravitational  force  on  the  drop, 
acoustic  levitation  is  achieved.  Several  features  of  acoustic 
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levitation  are  particularly  useful  in  this  experiment.  First  of  all, 
the  host  water  is  a  very  good  container  for  the  superheated  drop  as 
mentioned  above.  Secondly,  the  drop  can  be  held  in  water  very  steadily 
in  a  certain  position.  In  the  traditional  method  of  oscillating  drops, 
measurements  have  to  be  made  while  the  drops  are  rising  or  falling.  So 
acoustic  levitation  allows  accurate  measurements  of  the  drop's  size  and 
oscillation  frequency  with  comparative  ease.  Furthermore,  another 
acoustic  wave  can  be  coupled  into  the  system,  providing  a  means  to 
excite  the  drop  into  oscillations  without  touching  the  drop  with  any 
solid  object. 

The  acoustic  force  exerted  on  a  levitated  drop  is  a  known  function 
of  the  densities  and  compressibilities  of  the  drop  and  host  liquids.  By 
a  simple  technique  of  comparing  the  voltages  needed  to  levitate  two 
drops  of  different  liquids,  one  as  a  reference  and  the  other  one  as  an 
unknown,  at  the  same  position,  the  parameter  "compressadensity"  of  the 
unknown  liquid  can  be  obtained  providing  that  the  preperties  (density 
and  compressibility)  of  both  the  reference  and  host  liquids  are  known 
[5].  Compressadensity  is  a  known  function  of  density  and  adiabatic 
compressibility  as  expressed  in  Eq.  (60),  therefore,  if  we  know  the 
density,  then  the  compressibility  can  be  calculated.  This 
compressadensity  of  a  superheated  liquid  can  also  be  obtained  when  the 
interfacial  tension  is  measured.  It,  too,  is  difficult  to  measure  by 
other  methods  and  provides  additional  thermodynamic  data  in  temperature 
regimes  where  little  data  exist  (e.g.  [6]). 


1.5  Oscillations  of  a  Spherical  Drop 

To  evaluate  the  interfacial  tension  between  a  levitated  drop  and  its 
host  liquid,  we  need  a  model  relating  the  interfacial  tension  with  other 
parameters,  such  as  the  oscillation  frequency,  the  size  of  the  drop  and 
the  densities  and  viscosities  of  the  drop  and  host  liquids.  The 
equilibrium  shape  of  the  drop  should  be  approximately  spherical,  since 
the  gravitational  and  acoustical  forces  in  this  experiment  are 
insignificant  compared  to  the  interfacial  forces.  For  the  frequencies 
and  the  drop  deformations  in  this  work,  the  drops  behave  incompressibly . 
Therefore,  the  disturbance  of  the  interface  is  accompanied  by  the  drop's 
shape  oscillation  about  a  spherical  shape  without  volume  change.  The 
quadrupole  mode  (alternatively  prolate  and  oblate),  which  is  the 
simplest  mode  of  shape  oscillation,  is  studied  in  this  experiment. 

The  small  oscillations  of  a  liquid  mass  about  its  spherical 
equilibrium  shape  was  derived  by  Rayleigh  [7].  Lamb  generalized  the 
problem  by  supposing  that  the  liquid  drop  was  surrounded  by  an  infinite 
mass  of  another  liquid  [8] .  The  liquids  were  assumed  inviscid  by  both 
authors  so  that  viscosity  did  not  enter  their  expressions  for  the 
oscillation  frequency.  Lamb  also  discussed  the  effect  of  small 
viscosity  on  the  small  oscillations  of  a  mass  of  liquid  about  the 
spherical  form;  no  host  liquid  was  included  in  this  calculation  [8] . 
Miller  and  Scriven  [9]  did  an  analysis  of  small  oscillations  of  a 
viscous  fluid  droplet  immersed  in  another  viscous  fluid.  The  frequency 
and  rate  of  damping  of  free  oscillations  were  derived  for  arbitrary 
values  of  the  physical  properties  of  the  fluids  as  well  as  interfacial 
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viscosity  and  elasticity  coefficients.  Marston  [10]  independently 
derived  the  same  general  results  for  free  oscillations,  and  found  a  new 
term  to  correct  the  decay  time.  The  phase-frequency  relation  for  forced 
oscillation  was  also  established  for  the  first  time  in  this  work  of 
Marston. 

On  the  experimental  side,  Marston  and  Apfel  have  measured  the 
interfacial  tension  between  water  and  p-xylene  [11].  The  p-xylene  drops 
were  levitated  in  water  and  excited  into  quadrupole  shape  oscillations 
by  acoustic  forces.  The  oscillations  were  detected  by  a  rainbow 
photometry  technique  which  utilized  light  scattered  by  a  drop  at  the 
scattering  angle  normally  associated  with  the  rainbow.  The  interfacial 
tension  obtained  was  4%  lov/er  than  static  measurements  with  an 
uncertainty  of  about  10%.  Trinh,  Zwern,  and  Wang  [12]  have  done 
experimental  studies  of  drop  oscillations  with  emphasis  on  the 
frequencies  of  the  first  few  modes  and  the  damping  constant  for  the 
fundamental  (quadrupole)  mode.  They  used  a  slit  and  the  shadow  of  the 
uniformly  illuminated  drop  to  detect  the  oscillation. 

Viscosity,  which  causes  energy  dissipation,  should  relate  to  the 
damping  of  the  drop  oscillation.  In  his  theory,  Marston  gave  a  relation 
between  a  measured  damping  constant  and  the  viscosities  of  the  two 
liquids,  though  the  experiment  data  for  p-xylene  drops  in  water  did  not 
show  very  good  agreement  with  it  [11] . 


1.6  Fowkes1  Equation  for  Interfacial  Tension 

There  is  a  need  for  a  theory  that  predicts  the  interfacial  tension 
between  two  liquids  based  on  some  other  properties  such  as  the  surface 
tensions  of  each  liquid  against  its  vapor.  Unfortunately,  there  has  not 
been  a  theory  that  can  be  satisfactorily  applied  to  all  different  cases. 
Here,  we  briefly  introduce  a  semi-empirical  model  developed  by  Fowkes 
[13]  for  water-hydrocarbon  interfaces.  Fowkes1  model  has  been  widely 
used,  though  it  is  also  criticized  by  some  people  (e.g.  [14]). 


The  interfacial  tension  o.Q  between  liquid  A  and  B  is  the  work 
needed  to  expand  one  unit  area  of  the  interface  by  bringing  molecules  A 
and  B  from  the  interior  to  the  interface.  It  may  be  regarded  as  the  sum 
of  the  work  for  overcoming  the  cohesion  between  similar  molecules  of 
each  liquid  minus  the  work  done  by  adhesion  between  the  two  different 
species  across  the  interfacial  region.  So 


=  °A  +  °B 


where  W  is  the  work  done  per  unit  area  by  the  adhesion  or, 
equivalently,  the  increase  in  free  surface  energy  upon  separating  the 
interface  AB.  There  are  different  kinds  of  intermolecular  forces  such 
as  the  London  dispersion  force,  dipole-dipole  force,  and  hydrogen 
bonding.  For  the  water-hydrocarbon  system,  Fowkes  assumed  that  the 
interaction  between  water  and  hydrocarbon  molecules  is  only  due  to  the 
dispersion  force,  since  the  hydrocarbon  molecules  are  non-polar.  The 
geometric  mean  was  applied  to  give  W 


so  that 


where  the  subscripts  W  and  H  represent  water  and  hydrocarbon 
respectively,  and  the  superscript  d  stands  for  dispersion  part  of  the 

d  h 

intermolecular  force.  For  hydrocarbons  a„  =  o„.  For  water  ort  =  or,  + 

H  H  W  W 

ow^  where  the  superscript  h  denotes  the  intermolecular  force  due  to 

hydrogen  bond.  The  value  of  aj^  was  given  as  21.8±0.7  dyne/cm  according 

w 

to  Fowkes1  calculation  based  on  the  literature  data  of  the  interfacial 
tensions  between  water  and  each  of  eight  kinds  of  hydrocarbons. 

The  data  obtained  in  the  experiments  described  in  this  manuscript 
will  provide  additional  tests  of  the  applicability  of  different  theories 
of  interfacial  tension,  such  as  Fowkes'.  It  should  supplement  the 
existing  interfacial  tension  data  in  a  region  never  before  investigated, 
and  could  in  the  future  be  useful  in  the  development  of  more  accurate 
models  of  liquid-liquid  interfaces. 

1.7  Purpose  and  Scope 

The  purpose  of  this  study  is  to  obtain  experimentally  the  interfacial 
tension  between  water  and  each  of  three  superheated  liquids,  butane, 
isobutane,  and  propane,  at  room  temperature  and  atmospheric  pressure. 

The  data  and  method  of  this  work  may  have  applications  in  different 
areas.  For  example,  as  discussed  before,  they  may  be  important  in 
modeling  the  early  stage  of  the  phase  transition  (vapor  explosion)  of 
these  superheated  liquids  in  contact  with  water.  Also  the  data  may  be 
useful  in  testing  or  developing  models  on  interfacial  tension.  In 
addition  to  interfacial  tension,  measurements  of  the  compressibilities 


of  the  superheated  liquids  will  be  reported,  thereby  providing  data  on 
this  thermodynamic  property  previously  unavailable  in  the  literature. 
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The  technique  developed  for  this  experiment  is  based  on  the  shape 
oscillations  of  one  liquid  drop  in  another  liquid--the  host.  The  method 
developed  by  Marston  and  Apfel  [11]  has  been  modified  to  fit  this 
special  circumstance  and  to  improve  the  accuracy.  The  detection  of  the 
drop's  oscillation  is  adapted  from  the  work"  of  Trinh,  Zwern,  and  Wang 
[12]. 


A  simplified  model  for  the  quadrupole  oscillation  of  a  drop  in 
another  liquid  will  be  presented  in  chapter  2.  The  purpose  is  to  reveal 
the  physical  nature  of  the  problem  more  explicitly  while  appealing  to 
less  sophisticated  mathematics.  The  simplified  model  will  also  be 
compared  with  the  previous  theories. 

In  the  next  chapter,  we  begin  with  a  review  of  the  previous 
theories  on  the  shape  oscillation  of  a  drop.  A  simplified  model  is  then 
developed.  This  chapter  ends  with  the  equations  needed  for  the 
compressadensity  measurement  employing  the  established  acoustic 
levitation  technique.  In  chapter  3  and  4  we  describe  the  experimental 
equipment  and  procedures  respectively.  Then  we  present  and  discuss  the 
experimental  results  in  chapter  5.  Finally,  chapter  6  will  include  a 
summary  and  conclusions  of  this  work. 


Chapter  2 


THEORY 

In  this  chapter  we  discuss  the  theoretical  model  which  we  use,  in 
conjunction  with  experimental  measurements,  to  determine  interfacial 
tension.  We  first  describe  the  problem  and  basic  assumptions.  We  then 
review  several  previous  theories  for  the  oscillation  of  a  drop  in 
another  fluid.  Most  of  this  chapter  is  devoted  to  the  derivation  of  a 
simplified  model  which  we  developed  for  computing  the  damping  of  a 
freely  oscillating  drop  and  the  phase- frequency  relation  of  the  forced 
oscillation  of  a  drop  in  its  quadrupole  mode  (section  2.3).  The  chapter 
ends  with  a  brief  review  of  equations  used  in  the  acoustic  levitation 
technique  for  measuring  compressadensity . 

2.1  Problem  and  Basic  Assumptions 

The  physical  phenomena  of  the  method  of  oscillating  drops  for  measuring 
interfacial  tension  can  be  simplified  and  described  as  a  drop,  which  is 
immersed  in  another  fluid  and  has  a  spherical  equilibrium  shape, 
undergoes  shape  oscillation  with  small  amplitude,  under  the  following 
assumptions : 

1)  The  gravitational  force  and  the  levitational  acoustic  force  on 
the  drop  balance  each  other,  and  their  effects  on  the  equilibrium  shape 
of  the  drop  are  negligible  compared  with  the  interfacial  tension. 
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Therefore,  the  drop  has  no  steady  motion  relative  to  the  host  fluid,  and 
the  drop  is  spherical  at  equilibrium. 

2)  The  drop  behaves  incompressibly ;  thus  there  is  no  volume  change 
but  only  shape  change  accompanies  the  oscillation. 

3)  The  interface  is  free  from  surface  active  contaminants,  and  the 
interfacial  tension  is  a  constant,  uniform  value. 

4)  The  amplitude  of  the  shape  oscillation  is  so  small  that  the 
response  of  the  drop-host  system  is  linear  with  the  driving  fofce. 

5)  The  liquids  of  the  drop  and  the  host  are  immiscible,  so  that 
there  is  negligible  mass  transport  across  the  interface.  The  interface 
is  so  thin  (see  section  1.1)  that  it  may  be  considered  as  a 
two-dimensional  discontinuity;  the  physical  properties  of  the  fluids  on 
each  side  of  the  interface  are  uniform. 

2.2  Previous  Theories  for  Shape  Oscillations  of  a  Drop 

For  analyzing  the  above  problem  of  an  oscillating  drop  in  a  host, 
several  previous  theories  are  relevant  and  are  briefly  reviewed  in  this 
section.  We  first  review  some  theories  in  which  both  the  drop  and  host 
liquids  were  assumed  inviscid.  We  then  review  some  theories  in  which 
small  viscosity  of  both  the  liquids  were  assumed. 

2.2.1  Drop  and  Host  Liquid  of  Negligible  Viscosity 

2. 2. 1.1  Resonance  Frequency 
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Rayleigh  [7]  investigated  small  oscillations  of  a  liquid  mass  about  its 
spherical  shape  of  equilibrium,  and  obtained  the  resonance  frequency 
rn(n-l) (n+2)o 

«  =  I - 1  ,  (2) 

L  pR3  J 

where  u  is  the  angular  frequency,  n  is  the  mode,  a  is  the  surface 
tension,  p  is  the  density  of  the  liquid,  and  R  is  the  radius  of  the  drop 
at  equilibrium.  The  host  medium  and  viscosity  were  not  considered,  which 
may  be  an  acceptable  approximation  in  studying,  say,  a  water  drop  in 
air. 


Lamb  [8]  later  discussed  the  small  oscillations  of  a  liquid  drop 
about  its  spherical  form,  and  generalized  the  question  by  supposing  that 
the  drop  is  surrounded  by  an  infinite  mass  of  another  liquid.  He 
obtained  the  following  resonance  frequency: 
r  n(n+l)  (n-1)  (n+2)c  -yi 

u»*  =  |  -  |  ,  (3) 

L  {(n+l)p^  +  npQ}R3  J 

where  «*  is  the  angular  frequency,  n  is  the  mode,  o  is  the  interfacial 
tension,  p^  and  pQ  are  the  densities  of  the  inner  (drop)  and  outer 
(host)  liquid  respectively,  and  R  is  the  mean  radius  of  the  drop.  When 
pQ  is  negligible  compared  with  p^,  Eq.  (3)  is  equivalent  to  Eq.  (2). 

The  viscosities  of  the  two  liquids  were  ignored  in  the  derivation  of  Eq. 
(3).  The  calculated  values  of  w*  are  higher  than  the  measured  resonance 
frequency,  for  example,  by  about  10%  for  a  p-xylene  drop  of  mm  size  in 
water  [11],  where  both  liquids  are  similar  in  density  and  in  viscosity. 


2. 2. 1.2  Damping 


Lamb  [8]  also  discussed  the  effect  of  small  viscosity  on  the  small 
oscillation  of  a  liquid  drop.  Assuming  that  the  motion  is  irrotational 
and  the  velocity  potential  of  inviscid  flow  is  applicable,  he  calculated 
the  energy  dissipation  and  gave  the  damping  constant: 

1  R2 

-  =  -  ,  (4) 

s  (n-l)(2n+l)  v 

where  s  is  the  damping  constant  and  v  is  the  kinematic  viscosity  of  the 
liquid  of  the  drop.  The  damping  constant  is  defined  in  such  a  way  that 
if  A  is  the  amplitude  of  a  damped  oscillation,  t  is  the  time,  then  A  « 

e  tS.  Eq.  (4)  ignores  the  host  medium,  which  is  reasonable  if  the  host 
is  a  gas  of  negligible  density  and  viscosity. 

Lamb  also  derived  the  damping  for  a  spherical  bubble  in  a  liquid 
with  viscosity  v  by  a  similar  method: 

1  R2 

-  =  -  .  (5) 

s  (n+2)(2n+l)  v 

For  a  quadropole  mode,  n=2,  this  equation  gives  an  s  that  is  four  times 
the  s  of  Eq. (4) . 


2.2.2  Drop  and  Host  of  Viscous  Liquids 

Miller  and  Scriven  [9]  gave  a  rather  comprehensive  analysis  of  a  viscous 
drop  immersed  in  another  viscous  fluid,  by  which  the  frequency  and  rate 
of  damping  of  free  oscillations  can  be  calculated  for  arbitrary  values 


of  drop  size,  physical  properties  of  the  fluids,  and  interfacial 
viscosity  and  elasticity  coefficients. 


Marston  [10]  independently  derived  similar  results  for  free 
oscillation  of  a  drop  for  a  situation  in  which  both  drop  and  host 
liquids  have  small  viscosity  so  that  the  "viscous  wavelength"  at  each 
side  of  the  interface  is  much  smaller  than  the  radius  of  the  drop.  He 
found  a  new  term  to  correct  the  decay  time  for  free  oscillation.  Also, 
for  the  first  time,  Marston  derived  the  phase-frequency  relation  for 
forced  oscillation  of  a  drop  immersed  in  another  fluid.  In  his  theory, 
the  damping  constant  for  free  oscillation  is  giyen  as 

s  =  (41  )u*'2  +  -  ha2  .  (6) 

For  forced  oscillation,  the  phase- frequency  is  given  as 

tan£  =  +  yu)/(u*2  -  otu'^  -  u2)  .  (7) 

Here, 

(2n+l)2(vipopipo)ii 

0  =  — - - - - -  ,  (8) 

2'SRr[(piPi)i*  +  (u0p0)'5] 

(2n+i) [2(n2-l)yi2pi  +  2n(n+2)yQ2po  +  yiyQ[(n+2)pi  -  (n-l)pQ]} 

y  —  •  “■  ■  / 

(9) 

R2r[(yipi)4  +  <v0Hl2 

and 

r  =  (n+l)Pi  +  np0  .  (10) 

In  the  above  equations, 

s  =  damping  constant  of  the  free  oscillation  (see  Eq.(4)), 
u*  =  Lamb’s  natural  frequency,  as  given  in  Eq.  (3), 

5  =  the  phase  of  the  displacement  relative  to  the 

driving  force  in  the  forced  oscillation, 


w  =  the  frequency  associated  with  phase  £, 


n  =  the  mode , 


=  the  density  of  the  inner  (drop)  fluid, 

pQ  =  the  density  of  the  outer  (host)  fluid, 

=  the  dynamic  viscosity  of  the  inner  (drop)  fluid, 

yQ  =  the  dynamic  viscosity  of  the  outer  (host)  fluid,  and 

R  =  the  radius  of  the  drop  in  spherical  form. 

Marston  further  approximated  Eq.  (7)  by  applying  Newton's  method  of 
iteration  and  obtained  the  following  expressions: 

w(90° )  =  «*  -  (a/2)u*1"2  +  a2/4  (11) 

=  w*  -  s  , 

w(  5  )  =  w(90°)  -  s/ tan£  +  a2/(2tan£)2  (12) 

=  w(90°)  -  s/tan£  . 

Eqs.  (11)  and  (12)  were  used  directly  in  calculating  the  interfacial 
tension  in  the  experimental  work  done  by  Marston  and  Apfel  [11]  in  the 
following  way  (also  see  section  5.2).  The  damping  constant  s  was  given 
by  Eq.  (12)  with  measured  w(£)  and  w(90°),  and  w*  was  then  obtained  via 
Eq.  (11).  With  this  w*  substituted  in  Eq.  (3),  the  interfacial  tension  a 
was  obtained. 


2.3  A  Simplified  Model  for  Shape  Oscillation  of  a  Drop 

Eqs.  (2)  and  (3),  though  straightforward  in  derivation,  can  not 
accurately  describe  the  oscillation  of  a  liquid  drop  in  another  host 
liquid  due  to  the  neglect  of  viscosity.  On  the  other  hand,  Eqs.  (6)  and 
(7)  show  much  better  agreement  with  experiments,  but  the  derivation  is 
rather  lengthy  and  complicated.  With  the  purpose  of  revealing  the 


physical  nature  of  this  problem  more  explicitly  and  appealing  to  less 
sophisticated  mathematics,  we  developed  a  simplified  model.  The  basic 
idea  and  the  derivation  are  described  in  this  section. 

2.3.1  Analogy  with  the  Basic  Mass  and  Spring  System 

Although  the  motion  of  the  small  oscillation  of  a  drop  in  a  host  fluid 
in  a  certain  mode  is  three-dimensional,  there  is  only  one  degree  of 
freedom.  Let  us  compare  this  drop-host  system  with  a  basic  mass-spring 
system. 

2. 3. 1.1  Resonance  Frequency  of  a  Mass  and  Spring  System 

For  an  object  of  mass  m  acted  on  by  a  spring  of  spring  constant  k, 
without  driving  force,  the  equation  of  motion  is 
d2x 

m -  +  kx  =  0  , 

dt2 

where  the  x  is  the  displacement  of  the  mass  relative  to  its  equilibrium 
position.  The  natural  frequency  is  then 

u  =  (k/m)‘\  (13) 

When  this  mass-spring  system  is  subject  to  a  resistive  force, 
proportional  to  velocity,  and  is  driven  by  a  force  F0cos(u>t),  the 
equation  of  motion  becomes: 
d2x  dx 


m -  +  b -  +  kx  =  F0cosut. 

dt2  dt 


(14) 


The  resonance  frequency,  if  we  define  it  as  the  frequency  with  90°  phase 
angle,  has  the  same  form  as  that  given  in  Eq.  (13): 

w(90° )  =  (k/m)1*  .  (15) 

Notice  that  the  damping  term  b  does  not  appear  in  Eq.  (15),  therefore, 

the  spring  constant  k  can  be  obtained  by  knowing  only  w(90°)  and  m.  The 

phase-frequency  relation  may  be  illustrated  more  explicitly  by  an 

schematic  impedance  diagram  shown  in  Fig.  2.1(b),  accompanied  by  a 

vector  representation  of  Eq.  (14)  in  Fig.  2.1(a),  and  with  a  typical 

phase-frequency  diagram  in  Fig.  2.1(c)  [15].  In  Fig.  2.1(b),  the 

component  due  to  damping,  jwb,  has  a  phase  angle  of  90°;  this  fact  is 

responsible  for  the  resonance  frequency,  w(90°),  being  independent  of 

the  damping  coefficient  b. 


assume  x=A»cos(wt+£) 


(viscous  damping) 


Figure  2.1:  (a)  Geometrical  representation  of  Eq.  (14). 

(b)  Impedance  diagram  for  the  system  described  by  Eq.  (14). 

(c)  Phase  vs.  frequency  for  the  system  described  by  Eq.  (14) 
w(90J)  is  given  by  Eq.  (15). 

(d)  Schematic  impedance  diagram  for  the  shape  oscillation 
of  a  drop  in  another  liquid. 


2. 3. 1.2  The  Shape  Oscillation  of  a  Drop  in  Another  Liquid 


Let  us  consider  the  analogy  between  a  mass-spring  system  and  a  drop-host 
system.  The  interfacial  tension,  which  we  want  to  measure  in  the 
experiment,  is  equivalent  to  the  spring,  providing  the  restoring  force. 
The  equivalent  mass  of  the  drop-host  system  should  be  a  function  of  the 
densities  of  the  two  fluids  and  the  size  of  the  drop.  The  damping  is 
due  to  the  viscous  dissipation  in  the  drop  and  host  fluid. 

For  the  free  oscillation  of  a  drop-host  system  without  damping,  the 
natural  frequency  is  given  by  Eq.  (3).  For  the  quadrupole  mode,  n=2,  we 
have 


24o 


1% 


L  R3(3Pi+2po)  J 


(16) 


Eqs.  (16)  and  (13)  show  clearly  the  analogy  between  the  drop-host  system 
and  the  mass-spring  system.  However,  when  the  viscosity  of  the  tv/o 
fluids  is  non-negligible ,  due  to  the  nature  of  the  viscous  damping  in 
liquid,  the  resonance  frequency  w(90°)  is  no  longer  the  same  as  the  w* 
in  Eq.  (16).  The  u(90°)  is  not  only  a  function  of  stiffness  k  and  mass 
m,  but  also  a  function  of  viscosity, 

u(90°)  =  function  (k,  m,  viscosity). 


The  above  function  has  to  be  derived,  which  is  the  essential  goal  of  our 
simplified  model.  The  impedance  diagram  for  a  drop-host  system  may  be 
something  like  the  one  shown  in  Fig.  2.1(d). 


2.3.2  Velocity  Profile  for  an  Inviscid  Drop  in  an  Inviscid  Host 

For  simplicity,  all  the  following  expressions  are  for  oscillations  in 
the  quadrupole  mode  only.  Nevertheless,  the  method  can  be  applied 
similarly  to  other  modes.  Spherical  coordinates  with  polar  angle  0  and 
radius  r  are  used  to  describe  the  motion.  Due  to  the  symmetry  with 
respect  to  the  vertical  axis  through  the  center  of  the  spherical  drop, 
expressions  are  independent  of  the  azimuthal  angle. 

For  small  quadrupole  oscillation  about  a  spherical  form,  the  radial 
position  of  the  interface,  r,  at  any  instant  may  be  given  as  [8]  [7]; 

r  =  R„  +  eY2(8)sin(wt).  (17) 

where 

r  =  the  radial  position  of  the  interface  relative  to  the  center 
of  the  equilibrium  sphere, 

R  =  the  equilibrium  radius, 

R0  =  a  value  close  to  but  not  equal  to  the  equilibrium  radius  R 
r  4e2  -|-1 

=  r  |1  +  - j  ,  (18) 

L  5R2  J. 

t  =  the  amplitude  of  the  oscillation, 
u  =  the  angular  frequency, 

Y2(0)=  spherical  harmonic  function  of  order  2 
=  3cos20-1. 

The  spherical  harmonic  function  of  order  n  satisfies  the 
Laplace  equation  in  the  following  way: 

72(r\ )=0,  72(r'n'1Y  )=0. 


Assuming  the  flow  is  inviscid,  we  have  the  velocity  potential  which 
is  determined  by  the  boundary  condition  in  Eq.  (17),  at  the  internal 
points  [8], 

uR  r2 

=  -  —  e  Y2(0)  cos(ut)  ,  (19) 

2  R2 

and,  at  the  external  points, 
wR  R3 

*o  =  -  2  Y2(0)  cos(ut)  ,  (20) 

3  r3 

We  use  subscripts  i  and  o  to  denote  the  terms  related  to  the  inner 
(drop)  and  the  outer  (host)  flow,  respectively. 

The  corresponding  velocity  distribution  is  readily  obtained  from 
the  velocity  potential, 

u  =  ur®*r  +  =  -V0  ,  (21) 

where  the  subscripts  r  and  0  denote  the  components  in  the  direction  of  r 
and  0  respectively.  From  Eqs.  (19),  (20),  and  (21),  the  radial  and 
tangential  component  of  the  flow  velocity  in  the  inner  and  outer  region 


are : 

Uri  =  U/RWOcos^-DcosUt)  -  (22) 
u0i  =  (-3/ 2R) Eu>rsin(20) cos (ut)  ,  (23) 
Uro  =  (R/r)4Ew(3cos20-l )cos(ut)  ,  and  (24) 
u0Q  =  (R/r)4Ewsin(20)cos(wt)  .  (25) 


Fig.  2.2  shows  a  schematic  diagram  of  the  streamlines  at  a  certain 
instant  when  the  drop  is  moving  from  a  spherical  shape  toward  a  prolate 


For  very  small  amplitude  s,  the  velocity  at  the  interface  may  be 
considered  as  the  value  of  the  above  functions  evaluated  at  r=R.  Notice 
that  the  radial  velocity  is  continuous  across  the  interface  as  expected, 

ar-| 

—  I  =  uri  =  uro  =  eu(3cos2e-l)cos(u»t)  ,  (26) 

3t|r=R 

where  U  .  and  U  are  the  values  of  u  .  and  u  evaluated  at  r=R 
n  ro  n  ro 

respectively.  But  the  tangential  velocities  at  the  two  sides  of  the 

interface  are  opposite  in  direction  and  have  their  magnitudes  related 

as : 

uei  -  -  <3'2>"eo  (27) 

=  -  (3/2)EUsin(28)cos(wt)  , 

where  U  .  and  U  are  the  values  of  u,, .  and  u„  evaluated  at  r=R 
01  60  01  0o 


respectively. 
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2.3.3  A  Modified  Velocity  Distribution  with  Small  Viscosity 

The  velocity  distribution  predicted  by  the  velocity  potential  of 
inviscid  flow,  as  shown  in  Fig.  2.2,  is  found  to  agree  with  the  outcome 
of  a  qualitative  photographic  study  done  by  Trinh,  Zwern,  and  Wang  [12] . 
On  the  other  hand,  the  discontinuity  of  the  tangential  velocity  across 
the  interface  can  not  be  true  physically.  This  velocity  "jump"  has  to 
be  smoothed  out  by  viscosity. 

We  assume  the  viscosities  of  both  inner  and  outer  fluid  are  small, 
so  that  the  boundary  layer  thickness  (or  viscous  wavelength)  at  each 
side  of  the  interface  is  smaller  than,  say,  1/10  of  the  radius  of  the 
drop.  Then  we  may  model  the  flow  distribution  in  an  ad  hoc  way  by 
combining  a  thin  boundary  layer  at  each  side  of  the  interface  with  the 
potential  flow,  derived  on  the  inviscid  assumption,  at  all  points  other 
than  the  boundary  layer.  Based  on  this  approximate  flow  distribution, 
we  compute  the  damping  rate  and  the  phase-frequency  relation. 

To  consider  the  flow  in  the  boundary  layer  around  the  interface,  we 
take  a  rather  simple  approach.  We  first  estimate  the  tangential 

velocity  of  the  interface,  ug,  and  then  consider  the  boundary  layer  at 

each  side  of  the  interface  separately  as  a  consequence  of  an  enforced 
oscillation  of  a  solid  boundary. 

2. 3. 3.1  Tangential  Component  of  the  Interfacial  Velocity 


If  we  assume  that  there  is  no  pure  rotation  associated  with  this 
oscillation,  then  the  interface  should  oscillate  in  phase  with  the 


potential  flow,  as  does  u  --the  tangential  component  of  the  velocity  of 

•  s 

the  interface.  Physically,  we  know  that  the  amplitude  of  Ug  must  be  in 
between  U  .  and  U  .  The  ratio  U  /U  .  or  U  /U  should  be  determined  by 

8 1  qO  S  q1  S  00 

the  following  two  factors:  1)  The  momentum  carried  by  the  "external 
stream",  which  is  a  function  of  the  density  and  velocity.  2)  The 
diffusivity  of  momentum,  that  is,  the  viscosity. 

Assume  that  for  the  tangential  velocity  in  the  boundary  layer  the 
gradient  in  radial  direction  is  much  greater  than  the  gradient  in 
tangential  direction.  In  other  words,  the  boundary  layer  thickness  is 
much  smaller  than,  say,  a  quarter  of  the  circumference  of  the  space.  * 

Then,  locally,  we  may  approximate  the  flow  in  the  boundary  layer  as  a 
stream  of  viscous,  incompressible  fluid  over  a  parallel  stream  of 

1 

different  density  and  viscosity.  Also,  to  a  first  approximation,  we  may 

use  U„ .  and  U„  as  the  "external  stream  velocities".  Under  these 
01  0o 

approximation,  the  results  of  the  work  "  The  Velocity  Distribution  in 

m 

the  Laminar  Boundary  Layer  between  Parallel  Stream"  done  by  Lock  [16] 
can  be  applied  to  compute  the  Ug  as  follows: 

u.  X+( 6/y) 

c  =  -  =  -  ,  (28) 

u9i  1  + 

M 

6  r  1.425X  +  0.858  -,4 

-  =  |  -  |  .  (29) 

y  L  2 . 283Xpy  J 


A  V. 


“80  ‘  °o  Po  Pq 

where  X  =  -  =  -  -  (see  Eq.  (27)) ,  6=—,  p=— ,  y=—  , 
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6  denotes  boundary  layer  thickness,  p  denotes  density,  and  y  denotes 
dynamic  viscosity.  The  value  c  given  in  Eq.  (28)  depends  upon  the 
velocity,  density,  and  viscosity  of  the  two  parallel  streams,  as 


expected*.  With  a  given  Ug,  the  boundary  layer  at  each  side  of  the 

interface  can  be  modeled  as  the  flow  over  an  oscillating  boundary,  as 
described  in  the  next  subsection. 


2.3.4  Damping  of  Free  Oscillation 

In  this  subsection  we  discuss  the  damping  rate  of  the  free  oscillation 
of  a  drop  immersed  in  another  fluid,  making  use  of  the  flow  distribution 
described  in  section  2.3.3.  The  resultant  damping  rate  is  expressed  in 

terms  of  known  parameters  including  the  inner  and  outer  density,  Pi  and 

pQ  ,*  the  inner  and  outer  viscosity,  y^  and  yQ  ;  the  radius  of  the  drop 

R;  and  the  oscillation  frequency  w.  The  result  will  be  used  in  the 
derivation  of  phase-frequency  relation. 


*  We  give  two  extreme  cases  as  examples  for  Eq.  (28) 

a)  For  a  water  drop  in  air,  X=- (2/3 ) ,  pv=(p0v0)/(piPi)=2xlO-5, 

(6/y)>>l,  then  c=1.0,  which  means  that  U  is  very  close  to  U  . ,  the 

S  gl 

tangential  velocity  of  the  water  surface. 

b)  For  an  opposite  case  of  an  air  bubble  in  water,  X=-(2/3),  py=6xl04, 

(6/y)<<l.  c=X=-(2/3),  which  means  that  U  is  very  close  to  U.  ,  the 

s  00 

tangential  velocity  of  the  surface  of  the  water  host. 


The  damping  rate  is  a  function  of  the  instantaneous  total  energy  of 
the  system  and  the  rate  of  dissipation  at  that  instant,  as  governed  by 
the  following  expression: 
d 

—  (total  energy)  =  (rate  of  dissipation)  . 
dt 

In  the  problem  of  an  oscillatory  drop  in  a  host,  the  total  energy  T.E. 
is  the  sum  of  the  instantaneous  potential  energy  P.E.  and  kinetic  energy 
K.E.  The  potential  energy  is  a  function  of  the  interfacial  -tension,  and 
the  change  of  the  interfacial  area.  The  kinetic  energy  is  a  function  of 
the  density  and  the  flow  distribution.  The  energy  dissipation  is  due  to 
the  viscous  dissipation  in  the  boundary  layer  near  the  interface  and  in 
the  potential  flow  at  the  other  points. 

The  viscous  dissipation  may  be  considered  as  the  sum  of  the 
following  four  parts: 

1)  the  dissipation  in  the  inner  irrotational  flow, 

2)  the  dissipation  in  the  outer  irrotational  flow, 

3)  the  dissipation  in  the  inner  boundary  layer,  and 

4)  the  dissipation  in  the  outer  boundary  layer. 

2.3.4. 1  Dissipation  in  the  Inner  Irrotational  Flow 


The  dissipation  in  a  spherical  mass  of  fluid  with  radius  r,  calculated 
on  the  assumption  that  the  motion  is  irrotational,  is,  by  art.  355  of 
[8], 


(30) 
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2. 3. 4. 3  Dissipation  in  the  Inner  Boundary  Layer 

In  an  oscillatory  boundary  layer,  the  time  average,  over  one  cycle,  of 
the  rate  of  dissipation  in  the  boundary  layer  per  unit  area  of  the 
surface  (boundary)  is,  as  given  in  ['.!}, 

yU2 
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(34) 


where  U  is  the  amplitude  of  the  oscillatory  velocity  of  the  irrotational 
flow  outside  the  boundary  layer  (external  stream  velocity),  and  6  is  the 

boundary  layer  thickness.  The  boundary  layer  thickness  goes  up  with  the 

• 

viscosity  (diffusivity  of  momentum)  and  with  the  time  that  the  diffusion 
of  momentum  takes  place.  We  approximate  the  boundary  layer  thickness 
u 

as:  6=(2v/w)  [17].  Suppose  the  boundary  layer  thickness  is  much 

smaller  than  the  radius  of  the  drop,  so  Eq.  (34)  is  applicable  and  the 

amplitude  of  the  "external  stream  velocity"  may  be  approximated  by 

U„.-U  .  U  is  the  tangential  velocity  of  the  interface  and  U  =cU  .  (see 
0i  s  s  *  1  s  01 

Eq.  (28)).  On  making  use  of  Eqs.  (23),  (28),  and  (34),  integrati  .g  over 

the  interface  r=R,  we  have  the  rate  of  dissipation  in  the  inner  boundary 

layer  as: 


y(l-c)2[(3/2)£usin(20)]2 


12it(1-c)2u.  e2R2uT 


2(2v./w)  ‘ 


5(2V.)  ‘ 


2. 3. 4. 4  Dissipation  in  the  Outer  Boundary  Layer 


Again,  assuming  that  the  outer  boundary  layer  is  much  thinner  than  the 
radius  of  the  drop  and  the  amplitude  of  the  external  stream  velocity  may 

be  approximated  by  U  -U  ,  making  use  of  Eqs.  (25),  (27),  (28),  and 

00  s 

(34),  we  apply  the  same  method  as  that  in  the  above  paragraph  and  obtain 
the  time  averaged  rate  of  dissipation  in  the  outer  boundary  layer 


V  ( 1+1 . 5c) 2e2w2 (sin20 )2 


16it(1  +  1.5c)2ue2R2w~ 


2(2 Vj  «  ‘ 


15(2V_) 


2. 3. 4. 5  The  Total  Energy  of  the  Oscillation 


The  total  energy  at  any  instant  is  the  sum  of  the  instantaneous 
potential  and  kinetic  energy.  The  kinetic  energy  reaches  its  maximum 
value,  while  the  potential  energy  is  (set)  zero  (see  section  2. 3. 5. 2), 
at,  say,  t=0  in  Eq.  (17).  We  take  this  maximum  kinetic  energy  as  the 
total  energy  of  the  drop-host  system  in  oscillation.  Kinetic  energy  is 
a  function  of  velocity  and,  therefore,  is  a  function  of  velocity 
potential  [17], 

K.E.  =  hp  //0u«ndA 


=  hp  JI<t> urdA  . 

The  velocity  potential  with  an  assumed  amplitude  is  given  in  Eqs.  (19) 
and  (20).  Based  on  the  same  velocity  potential,  we  calculated  the  rate 
of  dissipation  in  the  previous  section.  Assuming  that  the  interface  is 
approximated  by  r=R,  and  substituting  Eqs.  (19)  and  (20)  in  the  above 
equation,  we  have 


K.E.)  .  =  (4/5)ire2R3p.w2 

inner, max 

Similarly,  we  have 


(37) 


K.E.)  .  =  (8/15)tte2R3p  w2 

outer, max  '  Ho 


Therefore,  the  total  energy  is 


T.E.)  =  (4/15)TTE2R3(3Pi  +  2po)w2 


(38) 


(39) 


2. 3. 4. 6  Damping  Constant 


If  the  total  energy  T.E.  =  ps2,  and  the  rate  of  dissipation  is  qe2,  then 


— (pe2)  =  -  qe2  . 
dt. 


This  reduces  to 
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r  -  r  e("4)(q/p)t 

t  "  £.j“  t 

where  e„  represents  the  amplitude  at  the  time  t=0,  and  e  represents  the 
amplitude  at  the  time  t.  If  we  define  damping  constant  S  in  such  a  way 

that  t  =  £0e  *  when  t=l/S<  then 

S  =  Os)(q/p)  •  (40) 

By  using  the  total  energy  and  the  rate  of  dissipation  obtained 
above,  together  with  the  Eq.  (40),  we  have  the  damping  constant  S  as: 


S  =  S,  +  S .  +  S  , 
b  l  o 


where 


(41) 


(2w)  2  i  j, 

Sb  =  -  [2.25(l-c)2(iiiP.)  2  +  (l  +  1.5c)2(vopo)  2], 

Rf 


15y. 


r2t 


40yo 


R2r 


r  =  3P.  +  2P  , 

Ki  Ko 


where  the  subscripts  b  denotes  the  contribution  of  dissipation  from  the 
boundary  layer,  i  denotes  the  contribution  from  the  inner  irrotational 
flow,  and  o  denotes  the  contribution  from  the  outer  irrotational  flow. 

Comparing  this  expression  with  Lamb's  formula,  we  see  that  the  is 


equivalent  to  the  s  in  Eq.  (4),  and  the  S  is  equivalent  to  the  s  in  Eq 
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2. 3. 4. 7  Examples  and  Comparison  with  the  Previous  Theories 

We  take  three  examples  to  illustrate  the  calculation  and  to  compare  the 
results  with  that  obtained  by  previous  theories  and  with  one  set  of 
experimental  data  from  previous  work. 

1)  For  an  air  bubble  of  radius  0.1  cm  in  water,  p ^=0.0012  g/cm3, 

y^=0 .018  cp,  po=0.997  g/cm3 ,  yQ=0.894  cp,  o=72  dyne/cm. 

(a)  By  Marston's  theory, (see  table  1  and  2  of  [10]) 

w*/2tt  =  148.03  Hz  (by  Eq.  (3))  , 
w  /2it  =  147.9  Hz  (by  Eq.  (11))  , 

V- 

s  =  Haw*  2  +  hZ  ~  Ha2 
=  0.6  +  17.75  -  0.0008 

=  18.35  Hz  . 

(b)  By  Lamb's  formula  (Eq.  (5)), 

v  =  Vq/Pq  =  0.897  cSt  , 
s  =  17.94  Hz  . 

(c)  By  this  simplified  model: 

w  /2-it  =  147.9  Hz  (a  given  value) 

(PoP0)/(PiV  =  4'lxl°4 

c  =  -  0.6667  (Us  =  U0o) 

S  =  Sv+  S.  +  S 
bio 

=  0.627  +  0.135  +  17.898 
=  18.66  Hz. 


The  damping  constants  obtained  by  both  models  are  very  close,  and  the 


simplified  model  gives  a  little  higher  value.  The  dominant  term  is  due 
to  the  dissipation  in  the  outer  irrotational  flow,  as  should  be 
expected. 

2)  For  a  water  drop  of  radius  0.1  cm  in  air, 

(a)  By  Marston's  theory, 

w*/2ir  =  120.87  Hz  (by  Eq.  (3)) 
u  /2ir  =  120.81  Hz  (by  Eq.  (11)) 

s  =  Haw*  2  +  H2f  “  Ha2 

=  0.376  +  4.498  -  0.00035 

=  4.87  Hz 

(b)  By  Lamb's  formula  (Eq.  (4)) 

v  =  =  0.897  cSt 

s  =  4.49  Hz 

(c)  By  this  simplified  model, 

w/2ir  =  120.81  Hz  (a  given  value) 

(Poyo)/(Piyi)  =  2-4x10’5 

c  =  1 .0  (u  =  u„ . ) 

s  01 

S  =  S,  +  S.  +  S 
bio 

=  0.378  +  4.480  +  0.241 
=  5.10  Hz  . 

Again,  the  calculated  damping  constant  obtained  by  this  simplified  model 
is  close  to  and  a  little  higher  than  ‘he  calculated  value  by  Marston's 
theory.  The  dominant  term  is  due  to  the  dissipation  in  the  inner 
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irrotational  flow,  as  one  expects. 

3)  Trinh,  Zwern,  and  Wang  did  some  experimental  work  on  the  free 
oscillation  of  a  drop  in  another  liquid  [12].  One  set  of  their  data  for 
phenetole  drops  in  a  mixture  of  water  and  methanol  was  compared  with  the 
theoretical  predictions  given  by  Marston's  theory  and  the  simplified 
model,  the  result  is  shown  in  Fig.  2.3.  It  shows  good  agreement  between 
the  data  and  the  calculated  values  based  on  both  models.  This  set  of 
data  is  chosen  because  of  the  small  viscosity  of  the  inner  liquid,  for 

which  R/6^  ~  10.  Good  agreement  is  also  found  between  other  sets  of 

data  and  the  simplified  model  even  when  R/fi^  ~  5,  although  the 

discrepancy  grows  as  the  ratio  R/6^  decreases,  as  we  expect.  For 
example,  in  the  experiments  with  drops  of  mixture  consisting  of  silicon 
oil  and  CCl4  in  water  host,  done  by  Trinh  et  al,  p ^=0.999  g/cm3, 

pQ=0 . 997  g/cm3,  p^=16.5  cp,  yQ=0.89  cp,  R-0.5  cm,  the  calculation  using 

our  simplified  model  gives  R/6^~3,  and  the  predicted  damping  constant  is 
about  40%  higher  than  their  measured  value. 


2.3.5  Phase-Frequency  Relation  in  Forced  Oscillation 


As  mentioned  in  section  2.2.1,  the  phase-frequency  relation  for  forced 
oscillation  is  readily  obtained  if  we  know  each  component  (both 
magnitude  and  phase)  of  the  impedance  corresponding  to  stiffness, 
inertia,  and  damping,  respectively.  Here,  the  impedance  is  defined  as: 

p  driving  pressure  (associated  with  quadrupole  mode) 

impedance  =  -  =  -  . 

A  displacement  (of  quadrupole  mode)  (42) 

Thus,  the  phase  of  the  displacement  is  relative  to  the  driving  pressure 

We  compute  each  component  of  the  impedance  separately  by  the 
following  approach.  First,  the  phase  angle  is  determined  by  physical 
reasoning.  We  then  compute,  by  energy  balance  or  force  balance,  the 
driving  pressure  associated  with  the  quadrupole  displacement  described 
by  r  =  R  +  eY2 (0)sin(ut)  (see  Eq.  (17)). 

2. 3. 5.1  Driving  Force 

In  experiments  described  in  this  manuscript,  the  drop  is  driven  into 
oscillation  by  a  modulated  acoustic  wave  (see  section  3.1).  Detailed 
analysis  of  the  acoustic  radiation  force,  which  is  a  time  averaged 
second  order  effect  [18],  exerted  on  the  interface  would  require  the 
first  order  pressure  distribution  of  the  incident  wave  in  the  specific 
configuration  of  this  experiment,  and  the  scattered  wave  incorporated 
with  the  boundary  condition  at  the  interface  of  the  drop  and  the  host. 
The  normal  and  tangential  component  of  the  resultant  ladiation  stress 
may  be  obtained  by  radial  and  tangential  projection;  further  expansion 


of  the  radiation  stress  in  terms  of  spherical  harmonic  functions  will 
give  the  normal  and  tangential  component  of  the  radiation  stress 
responsible  for  the  quadrupole  oscillation  [10]  [19]. 

Yosioka  and  Kawasima  [20]  derived  the  acoustical  radiation  force  on 
a  compressible  sphere,  the  tangential  radiation  stress  is  vanished  in 
their  derivation  due  to  the  inviscid  assumption.  Some  experiments  for 
measuring  "compressadensity"  using  acoustical  levitation  technique  is 
based  on  their  theory,  and  the  results  of  measurement  show  very  good 
agreement  with  the  theory  [5].  Herrey  [21]  has  shown  that  the 
tangential  radiation  stress  is  much  smaller  than  the  normal  stress  at  a 
highly  reflecting  metallic  surface  in  water  induced  by  a  sound  beam. 

For  a  liquid-liquid  interface,  the  boundary  is  relatively  "free"  to  move 
compared  with  the  boundary  of  liquid-solid  interface;  hence  the  velocity 
gradient  and  the  stress  at  the  liquid-liquid  boundary  would  be  expected 
even  smaller.  Also,  the  viscosities  of  both  inner  and  outer  liquids  in 
our  experiments  are  small.  In  the  following  calculation,  we  assume  that 
the  tangential  component  of  the  radiation  stress  is  much  smaller  than 
the  normal  component,  and  the  power  input  due  to  the  tangential 
radiation  stress  is  negligible  compared  with  the  contribution  of  normal 
radiation  stress. 

2. 3. 5. 2  Impedance  Due  to  Stiffness 

The  displacement  of  a  component  of  pure  stiffness  subject  to  a  driving 
force  is  in  phase  with  this  driving  force,  that  is,  £  =  0.  In  the 
problem  of  an  oscillating  drop,  stiffness  is  due  to  the  interfacial 
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tension,  as  assumed  in  section  2.1.  We  derive  the  impedance  due  to 
stiffness,  denoted  by  K,  in  two  different  ways  which  employ  force 
balance  and  energy  balance,  respectively,  as  follows: 


1)  The  fluid  pressure  is  discontinuous  at  a  curved  surface  of 
separation;  the  pressure  differential  is 


p.  -  Po  =  o(R1"1  +  Ra”1)  .  (43) 

This  is  Laplace's  formula,  where  Rt  and  R2  are  the  principal  radii. 

Since  we  are  only  interested  in  the  pressure  variation  accompanying  the 

quadrupole  oscillation,  let  us  set  pQ  =  0  at  equilibrium  (the  ambient 

pressure).  Then,  p^  =  o( 2/R)  at  equilibrium.  The  inner  pressure  p^ 

will  not  vary  with  the  shape  oscillation,  since  there  is  no  volume 
change  and  there  is  no  mass  involved  in  this  discussion.  The  interface 
is  described  as  (see  Eq.  (17)) 
r  =  R  +  EY2(0)sin(ut)  , 
which  gives  [8] 

1/Rl  +  1/R2  =  2/R  +  (4/R2)EY2(0)sin(wt)  .  (44) 

Therefore,  the  variation  of  pQ  accompanying  the  shape  oscillation  due  to 

the  interfacial  tension  is,  from  Eqs.  (43)  and  (44), 

Po  =  -  (4/R2)oeY2(0)sin(wt)  ,  (45) 

where  the  minus  sign  indicates  that  the  outer  pressure  at  the  two  poles 
reaches  a  maximum  when  the  drop  is  in  an  oblate  shape  with  maximum 
displacement.  For  this  reason,  we  define  the  displacement  A  in  the 
following  way,  so  that  the  displacement  is  in  phase  with  the  driving 
pressure  for  a  pure  stiffness. 


A  -  R-r  =  -  EY2(0)sin(ut)  • 


(46) 


Since  is  constant,  pQ  equals  the  driving  pressure  p.  Thus  the 

component  of  the  impedance  due  to  stiffness  is 
p  4o 

K  =  —  =  -  .  (47) 

A  R2 

2)  The  potential  energy,  set  zero  when  the  drop  is  in  a  spherical 
shape  at  equilibrium,  equals  the  product  of  the  interfacial  tension  and 
the  change  of  the  interfacial  area.  For  a  small  quadrupole  displacement 
with  an  amplitude  e,  as  in  Eq.  17  with  t=ir/(2u),  the  potential  energy  is 
(32-tr/5)oe2  [7].  The  driving  pressure  on  the  interface  responsible  for 
the  quadrupole  displacement  can  be  assumed  to  be  p  =  -PY2 ( 8) sin(ut ) , 
where  P  is  the  amplitude,  and  the  minus  sign  is  assigned  to  give  the 
correct  phase  angle  as  discussed  previously  (see  Eq.  (45)).  The 
instantaneous  power  input  per  unit  area  is  therefore 
dA 

p  —  =  PY22(8)Ewsin(wt)cos(u»t)  .  (48) 

dt 

The  integration  of  the  above  expression  over  the  interface  r=R,  and  over 
the  period  from  t=0  to  t=ir/(2u)  is  found  to  be  (8n/5)PeR2,  which  should 
equal  the  potential  energy,  (32u/5)oe2,  corresponding  to  the 
displacement  A  =  -eY2(9)  at  t=n/(2u>).  Hence  P  =  4oe/R2.  We  therefore 
have 

p  4o 
K  =  -  =  —  , 

A  R2 

which  is  identical  with  that  obtained  via  force  balance,  in  Eq.  (47). 


2. 3. 5. 3  Impedance  Due  to  Inertia 


The  acceleration  of  a  mass  is  in  phase  with  its  driving  force,  so  the 
displacement  lags  by  180°  the  driving  force;  £  =  180°.  Here  we  compute 
the  pressure  on  the  interface  required  to  drive  a  spherical  liquid  mass 
into  the  oscillation  described  as  r  =  R  +  eY2(0)sin(wt) .  The 
calculation  is  first  carried  via  force  balance  and  then  via  energy 
balance. 

1)  With  the  inviscid  assumption,  the  velocity  potential  is 
determined  by  the  boundary  condition,  as  given  in  Eqs.  (19)  and  (20). 
Then  the  pressure  variation  at  both  sides  of  the  interface  accompany 
this  potential  flow  is 
3*il 

Pi  =  pi  - 1  =  (1/2)piRw2EY2(8)sin(wt)  , 

at  | r=R 

9*ol 

po  =  pQ  - |  =  -(l/3)poR«2eY2(0)sin(wt)  . 

at  | r=R 

The  required  driving  pressure  is  p  =  p^  -  pQ,  and 

p  =  pi  -  p0  =  (l/6)Rw2(3Pi+2po)sY2(0)sin(wt)  .  (49) 

Thus  the  component  of  the  impedance  due  to  inertia,  denoted  by  M,  is, 
from  Eqs.  (42),  (46),  and  (49), 

P 

M  =  -  =  j2(l/6)w2R(3p^+2pQ)  , 


(50) 


It  is  worth  noticing  that  K  +  M  =  0  when  the  driving  force  has  a 
frequency  u>  that 
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2)  The  kinetic  energy  reaches  its  maximum  value 

(4ir/15)u2R3(3p^+2po)£2  (see  Eq.  (39)),  when  t  =  0,  and  it  is  zero  when 

t  =  n/(2u).  The  driving  pressure  may  be  written  as:  p  =  PY2(0)sin(wt) , 
so  that  the  acceleration,  d2A/dt2  =  £u>2Y2(0)sin(u>t) ,  is  in  phase  with 
the  driving  pressure.  For  instance,  when  t  =  ir/(2u),  the  driving 
pressure  is  positive  and  maximum  at  the  two  poles,  while  the 
acceleration  at  the  two  poles  are  inward  and  maximum.  Then,  during  this 
period,  from  t  =  0  to  t  =  n/(2w),  the  power  "input"  is  -(8h/5)PeR2  (see 
Eqs.  (46)  and  (48)).  This  power  input  should  equal  the  maximum  kinetic 

energy,  thus  -  (8tt/5)PeR2  =  -(4ir/15)w2ezR3(3p^+2po) .  We  obtain  an 

expression  for  P  and  therefore 
p  w2R 

M  =  —  =  j2  -  (3Pi  +  2Pq)  , 

A  6 

which  is  identical  with  that  obtained  via  force  balance  as  in  Eq.  (50). 


2. 3. 5. 4  Impedance  Due  to  Viscosity 

The  question  of  what  is  the  impedance  due  to  viscosity  is  equivalent  to 
the  question  of  what  is  the  driving  pressure  on  the  interface  needed  to 
supply  the  energy  dissipation  in  the  motion  r  =  R  +  eY2  ( 0 )  sin(u>t ) . 
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1)  Let  us  first  look  at  the  irrotational  flow.  In  the  irrotational 
flow  associated  with  the  oscillatory  drop,  the  velocity  at  every  point 
is  oscillating  in  phase,  since  the  fluids  behave  incompressibly .  The 
viscous  stress,  which  is  proportional  to  the  spatial  gradient  of 
velocity  distribution,  is  therefore  in  phase  with  the  potential 
velocity.  Hence,  the  driving  pressure  on  the  interface  is  in  phase  with 
the  potential  velocity  and  has  a  phase  lead  of  90°  over  the  displacement 
A  =  -eY2(8)sin(«t)  . 

The  radial  velocity  of  the  interface  is 
dA 

—  =  -euY2(6)cos(ut)  . 
dt 

The  driving  pressure,  which  is  in  phase  with  the  velocity,  is  assumed  in 
the  form 

p  =  -PY2(0)cos(ut)  • 

Then  the  instantaneous  power  input  by  the  driving  pressure  acting  on  the 
interface  is 
dA 

Up — dA  =  ( 16ir/5)  eujR2Pcos2  (wt)  ,  (52) 

dt 

which  should  equal  the  rate  of  energy  dissipation.  (a)  In  the  inner 
irrotational  flow,  the  rate  of  dissipation  has  been  obtained  (see  Eq. 
(32))  as: 

16Try^E2Ru2cos2  (wt)  . 

By  equating  this  rate  of  dissipation  with  the  power  input  (Eq.  (52)),  we 

get  an  expression  for  P,  P  =  5y.E/R,  which  leads  to  B  . ,  the  component 

x  pi 

of  impedance  due  to  viscous  dissipation  in  the  inner  irrotational  flow, 

P 

A  R 


B  . 
Pi 


(53) 


(b)  In  the  outer  irrotational  flow,  the  instantaneous  rate  of 
dissipation  is  (see  Eq.  (33)) 

(128/3)iryoE2Rw2cos2(wt)  . 

By  equating  the  power  input  (Eq.  (52))  and  this  rate  of  dissipation,  we 
have  B  ,  the  component  of  the  impedance  due  to  the  viscous  dissipation 
in  the  outer  irrotational  flow, 

Bpo  =  j(40/3R)you  .  (54) 

2)  Now  we  consider  the  energy  dissipation  only  in  the  boundary 
layer,  and  compute  the  associated  impedance.  The  velocity  distribution 
in  the  boundary  layer  may  be  described  as  [17] : 

v(y , t)  =  Uelut(l-exp[-(l+j)y/6] )  , 

where  Uelwt  is  the  external  stream  velocity,  y  is  the  coordinate  normal 
and  relative  to  the  boundary.  The  frictional  stress  at  the  boundary  is 

the  real  part  of  the  following  expression  [17]: 

au  |  U 

y —  |  =  y( 1+j ) -  eiwt  , 

ay  |y=o  6 

so  this  skin  friction  has  a  phase  lead  of  ti/4  over  the  external  stream 

velocity  Uelwt.  In  this  problem  of  an  oscillating  drop,  the  external 
stream  velocity  is  in  phase  with  the  potential  velocity,  as  discussed  in 

section  2. 3. 3.1.  The  tangential  component  of  the  acoustic  radiation 

stress,  which  is  the  driving  force,  though  small,  has  to  be  in  phase 

with  the  skin  friction,  in  considering  the  balance  of  tangential  stress 

on  the  interface.  Therefore  the  phase  of  the  radiation  stress  is  n/4 

ahead  of  the  potential  velocity  and  is  3n/4  ahead  of  the  interfacial 

displacement,  that  is,  £=135°.  For  the  motion  r=R+eY2 ( 0 ) sin(«t ) ,  we 
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have  (see  Eq.  (46)) 
dA 

—  =  -euY2(9)cos(«t)  . 
dt 

Writing 

p  =  -PY2(0)cos(ut  +  tt/4)  , 

the  power  input  averaged  over  one  cycle  along  the  whole  interface  is 
then 

2-ir/w  dA 

(w/2n)  S  dt  //  p—  dA  =  (4tt2^/5)eR2uP  ,  (55) 

t=0  dt 

where  the  contribution  of  the  power  input  due  to  the  tangential 
radiation  stress  has  been  neglected. 

In  considering  the  rate  of  dissipation,  we  assume  the  tangential 
radiation  stress  is  much  smaller  than  the  normal  component,  so  that  the 
model  of  the  flow  distribution  developed  in  section  2.3.3  is  still  a 
good  approximation  for  the  flow  in  the  forced  oscillation.  We  then  have 
the  time  averaged  rate  of  dissipation  in  thr  boundary  layer  as  given  in 
Eqs.  (35)  and  (36).  By  equating  this  time  averaged  rate  of  dissipation 
with  the  above  time  averaged  power  input  in  Eq.  (55),  we  have  an 
expression  for  P  and  thus  obtain  the  component  of  impedance  (see  Eq. 

(42))  due  to  viscous  dissipation  in  boundary  layer,  B^,  as: 

•3  /  ?  i,  u 

Bfa  =  [(j-l)//2](2/3)w^[2. 25(1-0)2^^)  2  +  ( 1  +  1 . 5c)  2  ( y^)  2]  .  (56) 

2. 3. 5. 5  Total  Impedance 

The  total  impedance  is  the  vector  sum  of  these  components  given  in  Eqs. 
(47),  (50),  (53),  (54),  and  (56). 


impedance  =  K  +  M  +  +  B, 


(57) 


where  B  =  B  .  +  B 

P  Pi  P° 


=  j(w/R)[5)ji  +  ( 40/ 3 )  ]  , 


K  =  4o/R2  , 


M  -  ( j2/6)u2R(3p^  +  2pQ)  , 

0/0  J,  \, 

Bfa  =  [(j-l)/^2](2/3)«  /^t2.25(l-c)2(y.p.) 2  +  (1  +1 . 5c) 2 (yopQ) 2] 

A  schematic  vector  diagram  of  the  four  components  is  shown  in  Fig.  2.4. 

We  obtain,  finally,  the  phase- frequency  relation  as  the  following 


(B_  +  B,  //2)  /  (K-M-B.//2) 


2. 3. 5. 6  Comparison  with  the  Theory  by  Marston 


Eq.  (7)  is  the  phase- frequency  relation  derived  in  the  work  done  by 
Marston.  We  change  the  form  of  Eq.  (7)  by  multiplying  each  term  with 

(3p^+2pQ)R/6  to  get  the  "equivalent  components  of  the  impedance"  as 

follows : 

K'  =  4o/R2  , 

M'  =  j2(l/6)w2R(3p.+2po)  ,  (59) 

j5w[yi2pi  +  (8/3)v02P0  +  (1/6)yiP0(4Pi"P0)] 

V  = - ; - 7 - ' 

R[(yiPi>”5  +  (y0p0)"2]2 

( j-1  )25u>3/2(iiiiuopipo)  2 

V  =  — - ; - 7-  • 

2  26[  (yipi) 2  +  (yQp0)  *1 

Compared  with  the  impedance  obtained  in  the  simplified  model,  they  have 
components  of  similar  dependence  upon  parameters  like  density, 
viscosity,  and  frequency  at  the  same  phase  angles.  The  components  due 
to  stiffness  and  inertia  are  identical.  For  both  inner  and  outer 

liquids  of  small  viscosity,  the  value  of  is  close  to  the  value  of 

B,  1 ,  whereas  the  value  of  B  is  much  greater  than  that  of  B  1 .  Two 
b  p  *  P 

examples  in  the  following  show  the  typical  numerical  values  of  impedance 

due  to  damping: 


a)  For  a  hexane  drop  in  water,  p^=0.655  g/cm3,  vtj_=0.29  cp,  pQ=0.997 

g/cm3,  yQ=0.89  cp.  On  making  use  of  Eqs.  (28)  and  (29),  we  have 

6/u=0.114,  c=-0.5,  therefore,  B  =0.15w3^2,  B  ,'=0.13w^2,  whereas 

o  o 

B  =0.13U/R,  B  *=0.07W/R. 

P  P 

b)  For  a  p-xylene  drop  in  water,  p^=0.857  g/cm3,  y^=0.614  cp.  Eqs.  (28) 

3/2 

and  (29)  give  6/p=0.19,  and  c=-0.4.  Therefore,  B^=0.22w 

B  ,'=0.17w3/2  ;  B  =0 . 15w/R,  and  B  '=0.05W/R. 
b  p  p 

Fig.  2.5  illustrates  the  phase- frequency  relation  for  a  hexane  drop  in 
water  predicted  by  Marston's  theory  and  the  simplified  model 
respectively  with  given  drop  size  and  properties  of  the  liquids 
(interfacial  tension,  density,  and  viscosity)  as  listed  in  Table  1. 

This  figure  is  meant  to  be  representative,  and  good  agreememt  between 
the  simplified  model  and  Marston's  theory  is  shown. 


Figure  2.5:  Phase  vs.  frequency  curves  given  by  the  simplified  model  and 
Marston's  theory  respectively:  Example  for  host,  water; 
drop,  hexane;  drop  radius,  0.0825  cm. 


2.4  Compressadensity  Measurement 


The  droplet-levitation  formula  is  given  by  [5] 
dP  2pg  B  *  p* 

P-  - - G( - -  -)  , 

dz  bs  es  P 

|1-P*/PI 

G  =  - : -  ,  (60) 

&s*/3s  '  (5P*-2p)/(2p*+p) 

where  P(Z)  is  the  acoustic  pressure  amplitude  of  the  standing  wave  along 

the  vertical  axis  (Z) ,  g  is  the  gravitational  constant,  $s*  and  are 

the  adiabatic  compressibility  of  droplet  and  host,  respectively,  and  p* 
and  p  are  the  densities  of  droplet  and  host  liquid,  respectively.  The 
nondimensional  function  G  is  called  the  "compressadensity"  function. 

The  above  expression  is  independent  of  the  size  of  the  drop,  based  on 
the  assumption  that  the  acoustic  wavelength  is  much  greater  than  the 
drop's  diameter.  On  the  other  hand,  adiabatic  compressibility  is  used 
when  the  thermal  wavelength  is  much  smaller  than  the  diameter  of  the 
droplet.  These  assumptions  hold  in  our  experiment  where  the  levitation 
frequency  is  about  52  kHz  and  the  radii  of  drops  are  smaller  than  1  mm. 

We  may  extract  information  about  a  liquid  by  using  a  method  of 
comparison  and  without  measuring  the  acoustic  field  [5].  A  drop  of  a 
"reference"  liquid  with  known  compressibility  and  density  is  first 
levitated  at  a  position  Z=Z0  with  a  levitation  voltage  V' ,  and  then  a 
drop  of  an  "unknown"  liquid  is  levitated  at  the  same  position  with  a 

voltage  V*,  The  pressure  amplitude  P(Z0)  and  pressure  gradient  (dP/dZ) 


are  assumed  each  proportional  to  the  levitation  voltage,  so  we  have  the 
result : 

P 1  (dP/ dZ) 1  rV'-,2  G(es'/Bs,  p'/p) 

-  =  l-l  = - -  (61) 

P*(dP/dZ)*  LV*J  G(3s*/es,  p*/p) 

where  the  superscript  1  denotes  the  reference  liquid,  and  the 
superscript  *  denotes  the  unknown  liquid.  Eqs.  (61)  and  (60)  give  us  a 
the  compressadensity  function  of  the  unknown  liquid  providing  that  p, 

B  p',  Ps'  are  known,  and  that  V'  as  well  as  V*  are  measured.  If  we 

know  p*,  then  3*  can  be  calculated  from  Eq.  (60),  (and  vice  versa). 


Chapter  3 
APPARATUS 

In  this  chapter,  we  focus  on  the  equipment  used  for  size  and 
phase-frequency  measurements.  A  schematic  diagram  of  the  equipment  is 
shown  in  Fig.  3.1.  We  briefly  mention  some  standard  instruments  used 
for  measuring  sound  speed  and  density  in  the  last  section. 

The  equipment  in  Fig.  3.1  may  be  described  in  four  parts.  The 
first  part  is  the  apparatus  used  to  levitate  a  drop  in  water  and  to 
drive  the  drop  into  shape  oscillations.  The  second  optically  detects 
the  drop's  shape  oscillations  and  converts  the  oscillations  into 
electric  signals.  These  signals  are  processed  by  the  third  part,  and 
the  frequencies  corresponding  to  several  phase  angles  are  measured.  The 
forth  part  takes  pictures  of  the  drop  and  some  calibrated  .reference 
objects.  These  pictures  are  used  to  measure  the  size  of  the  drop.  One 
section  of  the  following  will  be  devoted  to  describing  each  of  the  four 
parts  of  the  apparatus. 
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Figure  3.1:  Diagram  of  apparatus 
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3.1  Levitation  and  Shape  Oscillation 

The  levitation  cell,  which  holds  the  host  liquid  and  transducer,  is  a 
box  with  a  cross  section  of  5.4  cm  square  and  a  height  of  11  cm.  The 
box  is  made  of  Lucite  plates,  3.18  mm  in  thickness  for  the  walls  and 
6.35  mm  for  the  bottom,  and  it  is  open  at  the  top.  This  cell  sits  on  a 
lab  jack  which  adjusts  the  vertical  position.  The  transducer,  which 
converts  electric  energy  into  mechanical  motion  and  generates  waves  in 
the  host  water,  is  a  Channel  Industries  (Santa  Barbara,  CA)  5400,  Lead 
Zirconate  Titanate  piezoelectric  ceramic  disk  with  a  diameter  of  3.9  cm 
and  a  thickness  of  1.3  cm.  The  transducer  is  immersed  in  water  and  sits 
on  the  bottom  plate  of  the  cell. 

At  a  selected  height  of  the  water  in  the  cell,  and  with  an  input  to 
the  transducer  at  certain  frequencies  (e.g.,  two  frequencies  used  in  our 
work  are  about  52  kHz  and  510  kHz),  the  transducer-water-cell  system 
resonates  and  generates,  standing  waves  in  the  water  column  (see  section 
4.1).  A  drop  of  a  liquid  other  than  water  in  this  standing  wave  field 
is  subject  to  an  acoustic  radiation  force,  which  is  a  time  averaged  (dc) 
second  order  effect  of  the  wave.  This  force  tends  to  drive  the  immersed 
drop  toward  a  position  of  maximum  pressure  if  the  drop  is  more 
compressible  than  the  water  "host",  while  gravity  tends  to  push  the  drop 
up  or  to  pull  it  down  depending  on  the  densities  of  the  drop  and  the 
host.  If  the  drop  can  find  a  position  where  the  acoustic  radiation 
force  counterbalances  the  gravitational  force,  it  is  acoustically 
levitated  at  that  position.  The  radial  pressure  distribution  (maximum 
at  the  center)  forces  the  drop  to  the  center  of  the  water  column  if  the 


drop  is  more  compressible  than  the  host.  Also,  the  local  pressure 
distribution  affects  the  shape  of  the  drop,  while  the  interfacial 
tension  tends  to  keep  it  spherical. 

An  input  voltage  at  about  52  kHz  is  used  to  levitate  drops.  At 
this  frequency  and  a  proper  water  level,  drops  can  be  levitated  at  an 
input  of  only  a  few  volts  for  all  the  test  liquids  in  this  work.  At 
this  frequency,  the  wavelength  is  much  longer  than  the  diameter  of  a 
test  drop  (less  than  2  mm  mostly),  so  that  Eqs.  (60)  and  (61)  are 
applicable  for  calculating  compressadensity ,  and  the  equilibrium  shape 
of  tjie  drop  does  not  deviate  from  sphericity  much. 

At  about  510  kHz,  the  wavelength  is  comparable  with  the  diameter  of 
the  drop,  the  radiation  force  is  stronger  10],  and  the  deformation  of 
the  drop  is  greater.  We  can  further  enhance  the  deformation  by  turning 
this  signal  on  and  off  at  a  frequency  matching  a  resonance  frequency  of 
the  drop;  this  can  be  done  by  modulating  this  high  frequency  signal  (see 
Fig.  4.2). 

Therefore,  three  sinusoidal  electric  signals  are  combined  in  such  a 
way  to  give  v  driving  the  transducer,  that  is 

v^  =  v,  +  v  •  v 
t  1  c  m 

According  to  their  functions,  we  call  these  signals  the  levitation 
signal  (~52  kHz),  the  carrier  signal  (-510  kHz),  and  the  modulation 
signal  (~50  Hz),  associated  with  subscripts  1,  c,  and  m  respectively. 

The  circuits  for  multiplication,  using  a  MC  1595L  linear  four-quadrant 


multiplier,  and  for  summation,  using  a  ML  318S  integrated  circuit,  was 
built  by  Marston  in  previous  work  [11]  (see  Fig.  3.1). 

« 

The  levitation  signal  is  generated  by  a  versatile  circuit  built  by 
Baxter,  Apfel,  and  Marston  [22] .  This  circuit  has  features  that  are 
very  useful  in  tuning  the  acoustic  levitation  system.  We  can  manually 
scan  the  resonances  of  the  acoustic  levitation  system  and  view  the 
resonances  on  a  storage  oscilloscope  (Tektronix  Model  5103N) .  The 
ordinate  is  a  dc  manifestation  of  the  current  through  the  transducer, 
and  the  abscissa  is  linear  in  the  frequency.  A  Tektronix  P6016  current 
probe  senses  the  current  through  the  transducer  and  produces  a  voltage 
signal,  which  is  then  converted  by  the  versatile  circuit  to  the  dc 
manifestation. 

The  carrier  signal  is  generated  by  a  Krohn-Hite  Model  4300 
generator.  A  Krohn-Hite  Model  1000  generator  and  a  Kepco  Model  FG100A 
ramp  generator  are  used  for  the  modulation  signal.  The  Model  1000 
generator  operates  in  an  external  frequency  control  mode  and  is  driven 
by  the  ramp  generator,  thereby  providing  easy  manipulation  of  sweeping 
rate  and  direction. 

The  combined  signals  v  are  fed  into  a  Krohn-Hite  Model  7500 

amplifier  before  going  into  the  transducer.  After  amplification,  the 

rms  value  of  v^  is  about  4  V  or  less.  The  rms  value  of  v  is  about  5  V 

or  less.  A  Data  Precision  Model  5740  frequency  counter  (also  period 
meter)  and  a  Data  Precision  Model  3500  digital  voltmeter  are  used  in 
monitoring  the  system. 


3.2  Detection  of  Shape  Oscillation 

In  this  work,  we  study  small  amplitude  (in  the  order  of  10  ym)  shape 
oscillations  of  a  drop  (radius  in  the  order  of  1  mm)  which  can  not  be 

touched  with  any  solid  object.  These  two  features  suggest  the  use  of 

light  to  detect  the  shape  oscillation.  We  use  a  technique  adapted  from 

an  experimental  work  done  by  Trinh,  Zwern,  and  Wang  [12]. 

We  use  a  microscope  illuminator  as  a  light  source,  which  has  a 
Sylvania  1493,  6V,  20  W  bulb  with  a  short  filament  and  a  lens  with  a 
short  focal  length  to  give  an  approximately  parallel  beam.  A  dc  power 
supply  is  used  for  the  lamp.  The  light  beam  passes  horizontally  through 
the  levitation  cell,  then  is  enlarged  by  a  lens  and  is  projected  onto  a 
plate  with  a  slit  on  it.  The  slit  is  vertical  with  a  width  about  1  mm 
and  a  height  about  5  cm.  The  light  through  the  slit  is  focused  onto  the 
sensing  surface  of  a  photomultiplier  (ENI  type  9798B).  The  shadow  of  a 
levitated  drop  centers  across  the  slit  and  blocks  part  of  the  light 
going  through  the  slit.  When  the  drop  oscillates  in  its  quadrupole 
mode,  the  amount  of  light  going  through  the  slit  changes  accordingly, 
therefore,  an  ac  voltage  signal  coherent  with  the  drop's  oscillation  is 
obtained.  The  sensitivity  of  the  photomultiplier  is  a  function  of  the 
power  supply.  We  use  a  Pacific  Photometric  Instruments  Model  227  high 
voltage  power  supply  at  1000  to  1200  V  for  the  photomultiplier. 


The  shadow  of  the  drop  usually  is  not  a  uniformly  dark  circle,  but 
has  a  bright  area  in  the  center  part.  If  the  refraction  index  of  the 
drop  liquid  is  very  close  to  that  of  water,  as  the  index  of  refraction 
of  butane  is,  the  shadow  looks  like  a  ring  of  a  very  thin  band. 
Therefore,  we  use  a  adjustable  stop  to  block  the  light  inside  the  circle 
of  the  shadow.  Also,  we  use  an  iris  diaphragm,  which  is  concentric  with 
the  shadow  and  opened  a  little  bigger  than  the  shadow,  yielding  a  "halo" 
outside  the  shadow.  The  stop  and  the  iris  -diaphragm  reduce  the  light 
going  through  the  slit  and  result  in  two  advantages.  First,  the 
photomultiplier  can  operate  at  a  higher  sensitivity  to  give  a  stronger 
ac  signal  for  a  given  amplitude  of  a  drop's  oscillations.  Secondly,  the 
noise  due  to  light  scattering  by  particles  in  the  path  of  the  light  beam 
is  reduced. 

3.3  Frequency  and  Phase  Measurements 

We  need  to  measure  the  frequencies  corresponding  to  90°  and  some  other 
phase  angles  (of  the  drop's  oscillation  relative  to  its  driving  force) 
around  the  quadrupole  resonance  of  the  drop.  This  is  accomplished 
mainly  by  the  aid  of  a  lock-in  amplifier,  which  is  a  Princeton  Applied 
Research  Model  5101  lock-in  amplifier  with  a  common  mode  rejection 
typically  85  dB  at  1  kHz. 

Lock-in  amplifiers  have  an  output  channel  and  two  input  channels, 
one  channel  for  reference  and  the  other  channel  for  signals  to  be 
measured.  The  signal  is  compared  with  the  reference,  and  only  those 


signal  components  which  are  synchronous  with  the  reference  yield  a  net 
dc  output.  Noise  and  other  non-synchronous  signals  do  not  contribute  a 
net  dc  output,  but  only  ac  fluctuations  which  can  be  reduced  to  very 
small  value  by  time  averaging  (filtering).  There  is  a  phase  control  on 
a  lock-in  amplifier  which  allows  the  phase  of  the  reference  drive  to  the 
internal  circuit  of  phase  sensitive  detector  to  be  set  at  any  angle 
relative  to  the  input  reference  signal.  For  example,  a  sinusoidal 

signal,  say  v^cos (2nf ^t)  is  used  as  reference,  signals  to  be  measured 

have  a  synchronous  component  v^cos(2irfrt  +  £),  and  the  phase  control  is 

set  at  A-  Then  the  output  of  the  lock-in  amplifier  is  a  dc  voltage 
proportional  to  the  time  average  of  the  following  product: 

vrcos(2ufrt  +  A)  •  vicos(2nfrt  +  £)  . 

This  dc  output  reaches  maximum  when  A=£  ;  it  is  zero  when  A=£+90D.  It 
is  easier  to  take  measurement  at  zero-output  than  to  take  measurement  at 
maximum-output.  In  this  experiment,  by  setting  the  phase  dial  a=£+90°  , 
and  sweeping  the  modulation  frequency,  we  get  the  frequency  f(£)  which 
nulls  the  output  of  the  lock-in  amplifier.  This  process  of  sweeping  is 
viewed  on  the  storage  oscilloscope  with  the  output  of  the  lock-in 
amplifier  as  the  ordinate  and  the  sweeping  frequency  as  the  abscissa. 

We  use  the  Data  Precision  Model  5740  counter  to  measure  the  period  of 
the  modulation  signal,  then  take  the  reciprocal  of  the  period  to  get  the 
frequency,  because  it  is  easier  to  measure  period  than  frequency  in  the 
low  frequency  region  (on  the  order  of  100  Hz). 

When  the  Model  5101  lock-in  amplifier  operates  in  the  region  of 
about  100  Hz,  the  true  phase  shift  A  is  different  from  the  set  value  on 


the  dial.  So  we  calibrate  this  Model  5101  lock-in  amplifier.  The 
calibration  is  carried  out  by  applying  two  in-phase  signals  into  the 
reference  and  signal  channel,  followed  by  adjusting  the  phase  control 
for  zero  output.  Comparing  this  phase  reading  with  the  true  value 
(90°),  we  obtain  the  error  in  the  phase  reading  at  this  frequency.  The 
calibration  covers  the  frequency  range  from  40  Hz  to  140  Hz,  the  signal 
level  from  0.01  to  0.3  V,  and  concludes  that: 

the  phase  reading  from  the  dial  +  40°  =  the  true  phase 
with  error  less  than  0.2°  in  phase  angle.  Hereafter,  the  phase  set  on 
the  lock-in  amplifier  mentioned  in  this  manuscript  is  the  corrected 
phase. 

For  better  phase  accuracy,  the  amplitude  of  the  reference  signal  is 
set  to  1  V  rms  as  suggested  in  the  manual,  in  both  the  actual 
measurements  and  the  calibration. 

When  the  reference  signal  sweeps,  the  internal  reference  circuit  of 
the  lock-in  amplifier  automatically  tracks  but  with  some  delay.  A 
relation  between  the  slewing  rate  df/dt  and  the  phase  error  8  stated  in 

-3 

the  manual  is  df/dt  =  3«10  »f»0.  In  this  experiment,  f~100  Hz,  df/dt  £ 

1/30  Hz/sec,  so  the  error  in  phase  angle  introduced  by  sweeping  is  less 


3.4  Size  Measurement 


As  shown  in  Eq.  (3),  the  interfacial  tension  goes  with  the  cube  of  the 
drop's  radius.  For  a  drop  with  1  mm  radius,  an  acceptable  error  of  1% 
in  radius  will  cause  an  unacceptable  error  of  3%  in  interfacial  tension. 
The  size  measurement  was  the  most  significant  factor  of  the  experimental 
uncertainty  in  the  work  of  Marston  and  Apfel  [11]. 

In  this  experiment,  we  take  pictures  of  a  levitated  drop  and  some 
reference  objects,  then  enlarge  the  pictures  and  obtain  the  size  of  the 
drop.  We  use  a  Polaroid  #545  film  holder  loaded  with  a  Polaroid  Type 
55/Positive-Negative  4x5  land  film  positioned  right  in  front  of  the  slit 
and  facing  the  light  beam.  By  turning  the  toggle  switch  of  the  power 

supply  for  the  lamp  on  and  off  for  a  duration  about  half  second,  we  get 

an  exposure  of  the  shadow  of  the  drop.  The  shadow  is  about  ten  times 
larger  than  the  drop.  The  negative  is  then  further  enlarged  yielding  a 
positive  print  which  is  about  75  times  the  size  of  the  drop.  A  truss  is 
built  to  support  the  film  holder  firmly,  and  two  tracks  on  the  truss 
enable  the  film  holder  to  slide  in  and  out  the  optical  path. 

We  take  pictures  of  a  drop  before  and  after  the  frequency-phase 
measurement.  Then  we  take  pictures  of  the  following  reference  objects: 
three  precision  plug  gauges  of  diameters  1.395,  1.614,  and  2.018  mm 
(from  Zero-Check  Inc.).  When  we  take  the  picture,  the  gauge  is  located 
through  the  position  where  the  drop  was  levitated.  Care  is  exercised  to 

enlarge  these  negatives  by  the  same  factor.  The  average  of  the 

enlargement  ratios  (between  the  shadow  on  the  final  positive  print  and 


the  gauge  itself)  of  the  three  gauges  is  used  to  evaluate  the  drop's 
diameter.  The  consistency  of  the  three  ratios  indicates  the  accuracy  of 
the  size  measurement.  Generally,  they  are  good  to  within  0.15%. 

When  we  take  a  picture  of  each  gauge,  it  is  important  to  have  the 
gauge  located  through  the  position  where  the  drop  was  levitated.  A 
filar  telescope  with  magnification  of  20  and  a  ruler  horizontally 
affixed  to  the  plate  with  the  slit  on  it  serve  to  mark  the  position  of 
the  levitated  drop  in  three  dimensions. 

3.5  Measurement  of  Sound  Speed  and  Density 

We  use  standard  laboratory  equipment  to  measure  some  properties,  sound 
speed  and  density,  of  the  non-superheated  liquid  samples.  Since  these 
techniques  are  well  known,  we  will  not  discuss  these  methods  in  detail. 
The  values  obtained  are  incorporated  into  the  calculation  for 
interfacial  tension  or  compressibility. 

We  measure  the  sound  speed  in  hexane,  pentane,  and  heptane  by  using 
a  Nusonics  sonic  solution  monitor  Model  6105  (manufactured  by  Mapco 
Inc.,  Tulsa  Oklahoma)  (referred  to  Table  10).  The  container  for  the 
ultrasonic  probe  and  the  sample  liquid  is  placed  in  a  water  bath  which 
maintains  a  constant  temperature  with  variations  less  than  0.1°C.  The 
sound  speeds  obtained  using  this  technique  are  good  to  within  0.1%. 


For  measuring  density,  a  Mettler/Parr  DMA-40  densitometer  is  used. 
The  chamber  for  the  test  liquid  is  also  in  a  water  bath  to  maintain  a 
constant  temperature.  The  uncertainty  of  this  measurement  is  less  than 
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PROCEDURE 

The  procedure  for  measuring  interfacial  tension  in  this  experiment  is  as 
follows:  We  first  set  and  tune  the  levitation  system  to  get  a  strong 
acoustic  standing  wave  in  the  water  "host".  Then  we  put  a  drop  of  the 
liquid  to  be  measured  into  the  water.  While  the  drop  rises  up  due  to 
buoyancy,  the  acoustic  field  traps  the  drop  and  levitates  it  at  a 
certain  position.  We  take  pictures  of  the  enlarged  shadow  of  the  drop 
for  measuring  the  drop's  size.  Then  we  drive  the  drop  into  quadrupole 
shape  oscillation.  Around  the  resonance,  several  frequencies 
corresponding  to  different  phase  angles  of  the  oscillation  relative  to 
its  driving  force  are  measured.  Details  of  the  above  procedure  are 
described  in  the  following  sections. 

Some  other  parameters,  such  as  density  and  viscosity,  are  either 
measured  by  other  techniques  (see  section  3.5)  or  obtained  from  the 
literature.  The  procedure  for  measuring  compressadensity  by  the 
acoustic  levitation  technique  is  described  in  the  last  section. 
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4.1  Drop  Levitation 


The  host  liquid  in  the  levitation  cell  is  water  for  all  experiments  in 
this  study.  Surface  active  contaminants,  even  in  trace  amount,  may 
change  the  surface  or  interfacial  tension  significantly.  Therefore,  we 
are  very  careful  in  keeping  the  host  water  and  the  apparatus  components 
in  contact  with  the  host  water  clean.  Distilled  water  is  fed  into  a 
Barnstead  Nanopure-A  4-holder  system,  and  circulates  through  the  system 
to  reach  a  resistivity,  which  indicates  the  purity,  of  17.7  -  17.9 
megohm/ cm.  The  levitation  cell,  the  transducer,  and  others  like 
syringes  and  needles  are  kept  away  from  dirt  and  are  flushed  a  few  times 
with  the  clean  water  before  an  experiment  is  conducted.  The  clean  water 
is  degassed  and  then  slowly  poured  into  the  levitation  cell.  Then  we 
place  the  transducer  at  the  center  of  the  bottom  of  the  cell.  There  are 
no  bubbles  visible  in  the  cell  when  an  experiment  is  performed. 

We  tune  the  levitation  system  as  follows:  To  find  a  configuration 
at  which  the  transducer-water-cell  system  resonates,  we  adjust  the  water 
level  and  sweep  the  levitation  frequency,  while  monitoring  the  current 
going  through  the  transducer.  The  carrier  and  modulation  voltage  are 
zero  at  this  stage.  The  level  of  the  water  is  changed  by  taking  out  or 
injecting  in  small  quantities  of  water  with  a  syringe.  The  frequency  is 
swept  manually.  At  the  same  time,  we  view  the  resonances  of  the 
levitation  system  on  the  storage  oscilloscope  with  the  display  of 
current  versus  frequency  (see  section  3.1).  There  are  several 
configurations  of  water  level  and  frequency  which  give  sharp  peaks  on 
the  current-frequency  display.  Each  configuration  requires  a  different 


minimum  voltage  to  trap  a  drop  of  a  certain  kind  of  liquid,  and 
levitates  the  drop  at  a  different  position  with  this  minimum  trapping 
voltage.  The  water  level  (about  8.3  cm  from  the  top  of  the  bottom 
plate)  and  the  levitation  frequency  (about  52  kHz)  that  we  chose  provide 
stable  levitation  at  low  voltage  (about  4  V  for  hexane  drops).  The 
volume  of  water  in  the  cell  is  about  220  cm3.  The  levitation  position 
is  about  at  the  middle  of  the  height  of  the  water  column. 

To  levitate  a  drop,  we  put  it  into  the  bottom  part  of  the  cell 
(details  described  in  section  4.2)  and  trap  it  as  it  rises  (all  the 
liquids  to  be  measured  in  this  study  are  less  dense  than  water)  by 
turning  on  the  levitation  voltage.  After  tuning  the  levitation  system, 
we  levitate  a  hexane  drop  which  serves  for  checking  the  tuning  and  for 
aligning  the  optical  system.  A  stable  levitation  condition  at  about  4  V 
for  a  hexane  drop  indicates  good  tuning  of  this  levitation  system.  The 
positions  of  the  cell  and  the  enlarger  lens  are  adjusted  to  project  a 
sharp  shadow  onto  the  film  plane  and  to  center  across  the  slit.  Then  we 
use  the  filar  telescope  and  the  ruler  to  mark  the  position  of  the  drop 
in  three  dimensions  (see  section  3.4). 

After  finishing  these  procedure  for  preparation  ,  we  release  the 
hexane  drop  for  alignment  and  levitate  a  drop  of  the  liquid  to  be 
measured  at  the  same  position  as  that  of  the  hexane  drop  for  alignment. 
The  size  and  frequency-phase  measurements  are  then  conducted. 
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4.2  Drop  Making  and  Transportation 


For  non-superheated  liquids  like  hexane  and  pentane,  we  use  a  threaded 
plunger  syringe  and  a  needle  (from  Hamilton  company)  to  transport  the 
liquid.  A  drop  that  adheres  to  the  tip  of  the  needle  can  be  detached  by 
a  little  disturbance;  then  the  drop  rises  up  and  is  trapped  by  the 
acoustic  field. 


Introducing  superheated  drops  is  somewhat  complicated  compared  to 
non-superheated  drops.  The  process  is  shown  in  Fig.  4.1.  The  test 
liquid  is  contained  in  a  lecture  bottle  and  is  in  equilibrium  with  its 
vapor  at  room  temperature.  Some  of  the  liquid  is  drawn  into  a  threaded 
plunger  syringe  with  a  miniature  inert  valve  (from  Hamilton  Company) 
fitted  on  the  tip  of  the  syringe.  The  threaded  plunger  and  the  valve 
can  take  pressure  so  as  to  keep  the  liquid  in  equilibrium  (not 
superheated)  (referred  to  Fig.  4.1(a)).  The  syringe  is  precooled  for 
working  with  propane.  A  glass  tube  of  5  mm  i.d.  has  been  filled  with 
gel  and  cooled  in  a  freezer.  The  gel  is  composed  of  (by  weight)  1  part 
of  water,  1  parts  of  glycerine,  and  1.25  parts  of  water  soluble 
Aquasonic  100  ultrasound  transmission  gel  (from  Packer  Laboratories, 
Inc.,  NJ),  mixed  uniformly  and  centrifuged  to  remove  gas  bubbles.  We 
put  this  glass  tube  into  a  pressure  vessel,  pressurize  it  with  nitrogen 
gas  at  a  pressure  a  little  over  the  vapor  pressure  of  this  liquid  at 
room  temperature.  Then  we  push  the  needle  of  the  syringe  containing  the 
liquid  through  a  septum,  open  the  miniature  valve,  and  inject  the  liquid 
into  the  gel  forming  liquid  drops  (referred  Fig.  4.1(b)).  After 
releasing  the  pressure,  we  get  "superheated"  drops  in  the  gel.  The 


viscosity  of  the  gel  retards  the  drops's  moving,  and  the  smoothness  of 
the  gel  prevents  the  drops  from  vaporization.  This  technique  has  been 
adapted  from  a  technique  developed  by  Apfel  for  producing  neutron 
detectors  based  on  superheated  drops  in  a  gel  [23] . 

The  next  step  (referred  to  Fig.  4.1(c))  is  to  get  a  cluster  of  the 
gel  containing  one  superheated  drop  into  the  host  water  by  squeezing  a 
rubber  bulb  fitted  to  the  glass  tube.  This  cluster  of  gel  sinks  and 
sits  on  the  top  surface  of  the  transducer.  The  gel  gradually  dissolves 
in  water,  then  the  drop  floats  up  and  is  trapped  by  the  acoustic  field 
(referred  to  Fig.  4.1(d)).  The  amount  of  the  gel  going  into  water  with 
the  drop  is  no  more  than  0.1  g,  which  is  no  more  than  0.05%  of  the 
amount  of  host  water.  It  takes  about  20  minutes  for  the  gel  to  dissolve 
before  the  drop  is  released. 

All  the  apparatus,  such  as  the  syringe,  needle,  and  glass  tube, 
used  to  transport  the  liquid  are  cleaned  very  carefully  before  being 
used,  usually  flushed  by  acetone  followed  by  methanol  then  air  dried. 

We  always  use  the  liquid  directly  from  the  bottle  and  conduct 
measurements  immediately  to  reduce  the  possibility  of  contamination. 


Figure  4.1:  Procedure  for  introducing  superheated  liquid  drops  into  the 
host  water  in  the  levitation  cell. 


4.3  Procedure  for  the  Size  Measurements 

By  observing  the  drop  through  the  telescope  and  by  adjusting  the 
levitation  voltage,  we  can  align  the  drop  to  be  measured  at  the  same 
position  as  that  of  the  hexane  drop.  We  then  load  the  film  holder  and 
shoot  a  picture  of  the  drop's  shadow.  After  the  frequency-phase 
measurement  (section  4.4),  we  null  the  modulated  carrier  signal,  adjust 
the  levitation  voltage  to  move  the  drop  back  to  the  original  position, 
and  then  shoot  another  picture  of  the  shadow.  Then  we  release  the  drop, 
put  one  reference  gauge  into  the  cell  and  through  the  position  where  the 
drop  was  levitated,  and  take  a  picture  of  the  shadow  of  the  gauge.  Two 
other  reference  gauges  are  used  following  the  same  procedure.  It  is 
essential  for  accurate  size  measurement  to  position  the  reference  gauge 
through  the  levitation  position.  Very  good  stability  of  the  levitated 
drop  is  required  for  taking  good  pictures.  Therefore,  we  fine-tune  the 
levitation  frequency  before  exposures  of  the  drop.  These  negatives  are 
then  processed  for  measuring  the  size  of  the  drop,  as  mentioned  in 
section  3.4. 

4.4  Procedure  for  the  Phase-Frequency  Measurement 

After  taking  a  picture  of  the  levitated  drop,  we  slide  the  film  holder 
away  from  the  light  beam  and  adjust  the  stop  and  iris  diaphragm,  as 
mentioned  in  section  3.2.  The  light  going  through  the  slit  is  only  two 
bright  spots  immediately  outside  the  shadow  of  the  drop.  We  turn  on  the 
carrier  and  modulation  signal,  drive  the  drop  into  shape  oscillation, 
then  perform  the  frequency-phase  measurement. 


4.4.1  Tuning  of  the  Carrier  Frequency 

To  get  appreciable  amplitude  of  this  shape  oscillation,  the  carrier 
frequency  should  match  a  resonance  frequency  of  the  levitation  system  so 
that  enough  power  can  be  fed  into  the  transducer  to  generate  strong 
standing  waves  therefore  to  give  a  strong  radiation  force  on  the  drop; 
the  modulation  frequency  should  match  the  quadrupole  resonance  frequency 
of  the  drop  to  enhance  the  motion.  The  drop's  quadrupole  resonance 
frequency  can  be  estimated  from  the  apparent  size  of  the  drop  and  a 
priori  interfacial  tension.  The  carrier  frequency  can  be  tuned  by 
monitoring  the  response  of  the  drop  in  the  following  way: 

We  display  two  signals  on  a  Tektronix  T922  oscilloscope  (referred 

to  Fig.  3.1),  one  is  the  square  of  the  modulation  signal  (v^2)  as 

reference,  which  represents  the  driving  force  (explained  later  in  this 

subsection),  the  other  one  is  the  signal  (v  )  coming  from  the 

photomultiplier,  which  represents  the  drop's  response.  By  varying  the 
modulation  frequency  around  the  estimated  quadrupole  resonance  frequency 
of  the  drop  while  sweeping  the  carrier  frequency  (around  510  KHz),  and 

observing  the  two  traces  of  signals  (v^2  and  v  )  on  the  oscilloscope,  we 

can  get  a  signal  v^  that  appears  synchronous  with  the  reference  v^2 .  We 

fix  the  carrier  frequency  at  the  value  where  the  synchronous  signal 
appears  to  be  strongest  for  a  certain  modulation  frequency.  The  signal 

Vp  includes  noise  in  addition  to  this  synchronous  component.  A  low-pass 

filter  (Krohn-Hite  Model  3202  filter)  helps  us  visualize  this 


synchronous  component  when  the  output  of  the  photomultiplier  is  too 
noisy  (see  Fig.  3.1). 


We  confirm  by  the  following  observations  that  (the  synchronous 


component  of)  the  signal  v 


P 


corresponds  to  the  drop's  oscillation  in  its 


quadrupole  mode.  1)  The  signal  has  the  same  frequency  as  that  of  v^2 . 

2)  The  amplitude  of  the  signal  goes  up  or  down  in  phase  with  either  the 
carrier  or  the  modulation  voltage.  3)  The  phase  and  amplitude  of  the 
signal  vary  with  modulation  frequency  as  in  a  regular  resonance 
(referred  to  Fig.  2  1(c)).  4)  By  blocking  the  upper  or  lower  half  part 
of  the  light  beam  with  a  piece  of  card  board,  we  observe  that  the 
amplitude  of  the  signal  reduces  to  about  half,  and  the  phase  of  the 
signal  shows  no  change. 


The  frequency  of  the  forced  oscillation  of  a  drop  is  the  same  as 
that  of  the  exciting  force  (assuming  the  drop  is  a  linear  system),  but 
it  is  twice  the  modulation  frequency.  This  factor  of  2  comes  from  the 
fact  that  acoustic  radiation  pressure  is  a  time  averaged,  second-order 
effect.  A  schematic  diagram,  Fig.  4.2,  illustrates  this  frequency 

relation.  Therefore,  the  signal  vm2  with  an  amplitude  IV  in  rms  with 

respect  to  its  mean  is  fed  into  the  reference  channel  of  the  lock-in 
amplifier  (see  section  3.3  and  Fig.  3.1);  it  is  also  displayed  on  the 
Tektronix  T922  oscilloscope  to  help  in  tuning  the  carrier  and  modulation 

frequencies  (see  Fig.  3.1).  The  square  of  the  signal  v^  is  obtained  by 
feeding  v^  through  a  circuit  with  an  MC  1495L  linear,  four-quadrant 


multiplier. 


Figure  4.2:  Schematic  drawing  showing  the  time  variations  of  (a)  the 
carrier  signal  (~500  kHz),  (b)  the  modulation  s:.gnal  (~50 
Hz),  (c)  the  modulated  input  (v  x  v  )  responsible  for  the 
drop's  shape  oscillation,  (d)  acoustic  radiation  pressure, 
and  (e)  the  displacement  of  the  drop  with  a  phase  assumed  to 
be  90°. 
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4.4.2  Phase-Frequency  Measurement 

After  tuning  the  carrier  frequency  to  get  an  appreciable  amplitude  of 
the  drop's  oscillation,  we  feed  the  output  of  the  photomultiplier  into 

the  signal  channel  of  the  lock-in  amplifier,  apply  the  v^2  signal  to  the 

reference  channel,  and  set  the  phase  control  A  at  5+90°,  where  £  is  one 
of  several  selected  phase  angles.  The  dc  output  of  the  lock-in 
amplifier  is  displayed  on  the  storage  oscilloscope  as  the  ordinate, 
whereas  the  control  voltage  output  of  the  modulation  signal  (a  voltage 
proportional  to  its  frequency)  is  the  abscissa.  Then  we  sweep  the 
modulation  frequency  about  the  quadrupole  resonance  frequency  of  the 
drop,  and  obtain  the  f(£)  as  the  frequency  where  the  dc  output  of  the 
lock-in  amplifier  crosses  zero  (see  section  3.3).  We  sweep  the 
modulation  signal  very  slowly,  especially  when  the  trace  is  close  to  the 
zero-crossing  point.  Also  we  sweep  the  frequency  in  both  directions  to 
make  sure  that  the  sweeping  rate  is  slow  enough  that  sweeping  up  or  dov/n 
makes  negligible  difference  in  measured  f(£). 

We  repeat  the  process  for  several  different  phase  angles:  55°, 
67.5®,  90°,  112.5°,  and  125°,  which  are  centered  at  90°  and  chosen 
arbitrarily.  The  f(90°)  is  measured  at  least  three  times,  because  it  is 
the  most  significant  in  calculating  interfacial  tension.  A  drop  of 
either  superheated  or  non-superheated  liquid  can  be  acoustically 
levitated  for  more  than  two  hours.  But  a  levitated  drop  tends  to  be 
unstable  in  position  as  time  goes  on,  possibly  because  that  small 
variations  in  temperature  or  in  the  level  of  the  host  water  de-tune  the 
levitation  system  a  little,  or  because  small  particles  or  bubbles  are 


introduced  into  the  water  in  some  way.  Usually,  we  complete  the 
frequency-phase  measurement  in  30  minutes. 

4.5  Procedure  for  the  Compressadensity  Measurement 

The  measurement  of  compressadensity  may  be  either  incorporated  in  the 
interfacial  tension  measurement  or  as  a  separate  measurement.  We  first 
levitate  a  drop  of  liquid  A,  adjust  the  input  voltage  close  to  the 
minimum  trapping  voltage.  We  record  the  input  voltage  and  mark  the 
drop's  center  position  by  using  the  filar  telescope.  Then  we  release 
drop  A  and  levitate  a  drop  of  liquid  B.  We  vary  the  input  voltage  to 
bring  drop  B  to  the  same  position  as  was  drop  A,  recording  this  input 
voltage.  One  of  the  two  liquids  A  and  B  is  considered  the  "reference" 
with  known  density  and  compressibility.  By  applying  Eqs.  (60)  and  (61), 
we  can  calculate  the  compressadensity  of  the  other  "unknown"  liquid. 

The  host  liquid  is  water  for  all  tests.  The  diameters  of  the  drops, 
which  need  not  be  the  same  size  [5],  are  much  smaller  than  the  acoustic 
wavelength. 


Chapter  5 

RESULTS  AND  DISCUSSION 


In  this  chapter  we  present  the  results  of  the  quadrupole  oscillation 
experiments  described  in  chapter  4.  In  section  5.1  we  list  some 
relevant  physical  properties  of  two  common  liquids  (non-superheated)  for 
which  we  performed  experiments  in  order  to  test  the  accuracy  of  our 
measurements  for  interfacial  tension.  An  approximate  method  for 
reducing  data  used  in  a  previous  work  is  reviewed  in  section  5.2.  Some 
observations  and  problems  associated  with  experiments  involving  gel  are 
discussed  in  section  5.3.  A  method  for  inferring  inter facial  tension 
from  data  employing  a  least-squares  principle  is  presented  in  section 
5.4.  Section  5.5  deals  with  the  results  and  discussion  for 
non-superheated  liquids.  In  section  5.6  we  present  the  results  for 
superheated  liquids.  The  last  section  is  devoted  to  the  results  of  the 
compressadensity  measurements  for  three  superheated  liquids. 


5.1  Non-superheated  Sample  Liquids 

A  number  of  quadrupole  oscillation  experiments  were  performed  for  two 
common  liquids  in  order  to  test  the  accuracy  of  our  technique  for 
measuring  interfacial  tension  and  to  establish  the  method  for  data 
reduction.  N-hexane  and  n-pentane  were  chosen;  the  interfacial  tensions 
between  water  and  each  of  them  are  known  to  within  ±0.4%  (see  Table  1), 


and  they  are  saturated  hydrocarbon  compounds  similar  to  the  superheated 
liquids  studied  in  this  work. 

The  liquid  samples  were  from  Aldrich  Chemical  Company,  Inc.: 
n-hexane  with  99+%  purity,  catalog  #13938-6;  n-pentane  with  99+%  purity, 
catalog  #15495-4.  Some  relevant  properties  of  these  two  liquids  and 
water  are  listed  in  Table  1.  In  this  table,  the  densities  of  n-hexane 
and  n-pentane  were  measured  using  the  Mettler/Parr  density  meter,  and 
the  measured  values  agree  with  literature  values  (e.g.  [25])  very  well. 
The  viscosity  was  taken  from  literatures  [25],  as  was  the  density  of 
water  [26].  The  interfacial  tension  was  taken  from  the  literature  [27] 
[28]  and  is  compared  with  our  measurements.  In  this  article,  wherever 
hexane  and  pentane  are  mentioned,  they  are  abbreviations  for  n-hexane 


|  Liquid 

j  Temperature  | 

Density  | 

|  Viscosity 

|  Interfacial  Tension | 

1 

t  1 

1 

1  xlO2 

with  water 

deg  C 


n-hexane  j 

(CH3(CH2)„CH3)| 

- h 

n-pentane  | 

(CH3(CH2)3CH3)| 


0.6262 

0.6213 


0.240 

0.230 


dyne /cm 


51.1  ±  0.2 


1 - 1 

j  20  ! 

1 - 

|  0.9982 

-i — 
1 

1.002 

—1 - 1 

1  1 

|  25  | 

1  1 

0-9971 

» 

1 

i 

0.955 

1  1 
1  1 

Note:  for  water-hexane  interface,  do/ dt  =  -0.026  dyne/cm  #C,  where  o  denotes 
interfacial  tension  and  t  denotes  temperature  [24] . 


5.2  Approximation  by  Newton's  Iteration 
5.2.1  Approximation  Method 


Our  experiments  for  measuring  interfacial  tension  may  be  described 
briefly  as  follows:  We  drive  a  levitated  drop  into  forced  oscillation  at 
frequencies  around  its  quadrupole  resonance  and  measure  the  frequencies 
associated  with  several  phase  angles;  we  then  compute  the  interfacial 
tension  based  on  the  resonance  properties  or,  more  specifically,  the 
phase-frequency  measurements.  This  phase-frequency  relation 
incorporated  with  interfacial  tension  may  be  expressed  by  the  following 
equation,  according  to  the  model  derived  by  Marston  (see  Eq.  (7)), 

T  =  (oU3/2  +  Zu)/(ut*z  -  aw3/2  -  w2)  ,  (62) 

where  «*  is  the  inviscid  resonance  frequency  and  is  related  directly  to 

the  interfacial  tension  to  be  measured  as  in  Eq.  (16),  T  =  tan£,  £  is 

the  phase  angle  of  the  drop's  displacement  with  respect  to  its  driving 

force,  u  is  the  frequency  of  the  forced  oscillation,  a  and  y  are 

functions  of  the  drop's  size  and  of  the  properties  (density  and 

viscosity)  of  the  host  as  well  as  the  drop  liquid  (see  Eqs.  (8),  (9), 

and  (10)).  The  phase- frequency  relation  given  by  the  simplified  model 

derived  in  chapter  2  can  also  be  expressed  in  the  same  form  as  Eq.  (62), 

though  a  and  y  will  then  represent  different  functions  of  the  size  and 

the  properties  from  that  given  by  Eqs.  (3)  and  (9)  (see  Eq.  (58)). 

In  the  work  done  by  Marston  and  Apfel  [11]  Eq.  (62)  was 
approximated  in  the  following  way  and  then  employed  with  data  for 
deducing  interfacial  tension.  First,  Eq  (62)  may  be  rewritten  as  a 
function  of  u, 

f(w)  =  T«2  +  o(l+T)w3/^2  +  Yu  “  Tw*2  =  0  .  (63) 


By  applying  Newton's  method  of  iteration,  and  using  w*  as  an  initial 
approximation,  u  can  be  expressed  as  a  function  of  T  [24]  [10], 


f(w*) 

w  =  w*  -  -  (64) 

f'(w*) 

-  (w*  “  Haw*^  +  3/8a2)  -  1/T  (Haw* 2  “  3/4a2  +  hZ)  . 
Therefore, 


w(90° )  =  w*  -  Haw* 2  +  (3/8)a2  , 

(65) 

and 

AV  =  «(90«)  -  1/T  [Haw*1"2  -  (3/4)a2  +  HZ]  . 

(66) 

If  we  define  s  as 

s  =  [u($)  -  w(90° ) ]  (-T)  , 

(67) 

then  from  Eq.  (66),  we  have 

s  =  Haw*^  “  (3/4)a2  +  \Z  •  (68) 

The  damping  constant  s  can  be  obtained  experimentally  by  substituting 

the  data  of  phase-frequency  measurements  into  Eq.  (67),  and  sg  is  used 

to  denote  it  in  this  article.  For  an  interface  composed  of  two  liquids 
with  known  density,  viscosity,  and  interfacial  tension,  a  (from  Eq. 

(8)),  y  (from  Eq.  (9)),  and  w*  (from  Eq.  (16))  can  be  calculated,  and  a 

theoretical  prediction  of  the  damping  constant,  denoted  by  s^,  can  be 

obtained  via  Eq.  (68).  From  Eqs.  (65)  and  (68)  we  have 

w*  -  w(90°)  +  Haw***  -  (3/8)a2  , 

=  w(90° )  +  s  +  (3/8)a2  -  HZ  .  (69) 

Furthermore,  by  neglecting  the  last  two  terms  of  Eq.  (69),  we  have 

w*  =  w(90° )  +  s  .  (70) 

Eq.  (70)  was  used  to  infer  interfacial  tensions  from  measurements  in  the 

work  done  by  Marston  and  Apfel  [11].  In  this  way,  a  and  y  were  not 

involved  in  the  calculation  for  interfacial  tension. 


By  using  artificial  data  (that  is,  a  phase-frequency  relation  given 
by  Eq.  (62)  and  known  radius  as  well  as  properties  of  liquids),  we  can 
show  that  Eq.  (69)  is  a  very  good  approximation  for  calculating 
interfacial  tension,  and  Eq.  (70)  gives  higher  u*  than  Eq.  (69)  does  by 
about  0.7%  which  introduces  about  +1.5%  error  in  the  interfacial 
tension.  For  our  measurements  with  hexane  or  pentane  drops  in  water 
(those  marked  with  *  in  Table  5),  the  theoretical  values  given  by  Eqs. 
(8)  and  (9)  are  probably  good  approximations  for  a  and  y,  judging  from 

the  fact  that  the  measured  damping  constant  S£  is  very  close  to  the 

theoretical  value  s^.  Therefore,  we  estimate  that  the  error  in  deduced 

interfacial  tension  introduced  by  using  the  approximate  model  as  in  Eq. 
(70)  is  about  1.5%  too  high. 

5.2.2  Examples 

As  an  example  of  this  procedure  of  calculation,  we  list  the  data  for  one 
hexane  drop  in  Table  2.  To  deduce  the  interfacial  tension  from  these 
data,  we  first  apply  Eq.  (70)  and  have 
w*  a  (113. 69+4. 80)2n  =  744.50  , 
then  from  Eq.  (16), 

Oj  =  w*2R3I724  =  51.5  dyne/cm  , 

where  Oj  is  used  to  represent  interfacial  tension  obtained  by  using  Eq. 
(70).  With  the  literature  values  listed  in  Table  1,  the  theory  by 

Marston  gives  3^=27.64  (Eq.  (63)),  ctth=1.76  (Eq*  (8)),  and  2^=13.97 
(Eq.  (9)).  The  experimentally  determined  damping  constant  sg  (30.14)  is 
close  to  the  theoretical  prediction  s  ^  (27.64).  If  we  us  Eq.  (69)  with 


the  data  in  Table  2  and  the  above  ct^  and  y^,  we  get  the  corrected 

interfacial  tension,  denoted  by  o2 ,  as  50.6  dyne/cm.  This  set  of  data 
is  also  shown  in  Fig.  5.1.  In  this  figure,  triangle  circular  symbols 
represents  the  data  points,  and  the  dashed  line  represents  the 
theoretical  prediction  based  on  Marston's  theory  (Eqs.  (62),  (8),  (9), 
and  (16))  with  physical  properties  of  hexane  and  water  given  in  Table  1 
as  well  as  the  known  radius  of  the  drop.  The  data  points  follow  the 
slope  of  the  theoretical  line  closely,  which  indicates  agreement  between 

sg  and  s^.  The  measured  frequencies  associated  with  90°  and  other 

phase  angles  are  less  than  the  theoretical  prediction,  which  indicates 
that  the  measured  interfacial  tension  is  less  than  the  literature  value 
used  for  the  theoretical  prediction.  The  meaning  of  the  solid  line  will 
be  discussed  in  section  5.4.2.  This  is  a  typical  example  of  our 
measurements;  more  data  and  discussion  are  presented  in  section  ^.5  (see 
Table  5). 


TABLE  2 


Figure  5.1:  Phase  vs.  frequency  for  one  experiment  which  involved 
no  gel  in  the  preparation  of  the  drop. 


5.3  Gel  Effect 


As  described  in  section  4.2,  our  superheated  drops  was  first  formed  in 
gel  and  then  transported  into  the  water  host.  It  is  a  convenient  method 
for  introducing  superheated  drops  into  host  water;  other  possible 
alternative  methods  would  require  pressurizing  the  whole  levitation 
cell--up  to  10  atm  for  propane  drops  or  lower  pressure  with  some  cooling 
mechanism;  technically,  it  is  much  more  complicated  and  difficult. 
However,  because  of  the  contact  of  the  test  drop  with  the  gel  in  the 
process  of  transport  and  the  dissolution  of  gel  in  the  host  water,  one 
must  ask  whether  the  drop's  surface  is  contaminated  and  whether  the 
interfacial  tension  is  changed.  We  have  done  several  tests  to  explore 
this  issue. 

First,  we  measured  the  surface  tension  of  water  against  air,  then 
put  some  gel  (about  0.2%  by  volume)  into  the  water  and  measured  the 
surface  tension  again  after  the  gel  was  dissolved  in  the  water.  A 
Fisher  Surface  Tensiometer,  Model  20  (from  Fisher  Scientific  Company), 
employing  the  ring  detachment  method,  was  used  in  thf  measurements. 

There  was  no  detectable  difference  in  the  measured  surface  tension  (72.1 
±  0.1%  dyne/cm)  due  to  the  introduction  of  gel.  Secondly,  we  performed 
several  quadrupole  oscillation  experiments  for  hexane  and  pentane  (those 
marked  with  **  in  Table  5);  the  host  water  had  gel  dissolved  in  it 
(about  0.1  %  in  volume);  the  test  drop  was  injected  into  water  directly 
from  a  syringe.  The  measured  interfacial  tension  agrees  well  with  other 
data  which  were  taken  with  clean  host  water,  and  the  measured  damping 
constant  agrees  well  with  the  theoretical  prediction  (see  Table  5). 
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Therefore,  it  seems  that  the  gel  is  not  surface  active  in  water-hexane 
and  water-pentane  interfaces:  that  is,  there  is  no  tendency  for  the  gel 
to  pack  into  or  to  migrate  away  from  the  interface  and  to  change  the 
interfacial  tension  by  any  detectable  amount  [29] . 

On  the  other  hand,  for  the  experiments  of  hexane  and  pentane  which 
had  the  test  drop  first  formed  in  gel  (following  the  same  procedure  as 
in  Fig.  4.1  but  without  pressurization),  the  oscillation  of  the  drop 
exhibited  much  greater  damping;  the  measured  frequencies  at  selected 
phase  angles  showed  significant  difference  with  the  theoretical 
prediction:  the  frequencies  associated  with  90°  phase,  w(90°),  were 

lower  than  theoretical  prediction  by  about  5%,  and  the  values  of  s  were 

e 

higher  than  the  theoretical  value  s^  by  a  factor  about  2.  The  results 

of  these  experiments  are  marked  with  ***  in  Table  5.  One  example  is 
shown  in  Fig.  5.2;  the  circles  denote  data  points,  and  the  dashed  line 
represents  theoretical  prediction  as  mentioned  in  the  end  of  section 
5.2.  The  solid  line  will  be  discussed  in  the  next  section.'  These 
observations  suggest  that  there  might  be  some  residue  of  gel  on  the 
drop's  surface  while  the  measurements  were  taken,  though  it  could  not  be 
seen  from  the  pictures  of  the  drops.  These  observations  also  raise 
several  questions  about  how  to  analyze  the  measurements  which  have  gel 
involved:  First  of  all,  is  Eq.  (62)  still  a  good  model  for  the  motion? 
If  it  is,  then  what  is  an  appropriate  way  to  obtain  interfacial  tension 
from  the  measurements  of  resonance  properties?  The  values  of  a  and  JT 
need  to  be  estimated  in  order  to  employ  Eq.  (69)  or  to  assess  the 
validity  of  Eq.  (70).  The  fact  that  the  values  of  w(90°)  were  lower  and 


the  values  of  sfi  were  higher  than  theoretical  values  indicates  that  a 

might  be  greater  than  and  Z  might  be  different  from  the  theoretical 
values  given  by  Eqs.  (8)  and  (9)  (see  Eq.  (62)).  Therefore,  we  need  a 
new  scheme  to  analyze  data.  We  have  developed  a  method  employing 
least-squares  principle  to  estimate  a  ,Z,  and  interfacial  tension. 
Details  are  discussed  in  sections  5.4  and  5.5  and  the  Appendix  A  and  B 


5.4  Data  Reduction  by  the  Method  of  Least-squares 
5.4.1  Model  and  Procedure 


92 


Several  methods  have  been  considered  for  computing  interfacial  tension 
from  data  of  a  drop's  quadrupole  resonance  properties  (see  Appendix  A); 
the  one  we  chose  is  presented  here.  The  basic  idea  is  that  we  assume 
Eq.  (62)  is  an  appropriate  model  for  the  motion  but  leave  o,  X,  and  w* 
all  as  unknown  and  then  fit  this  model  to  data.  We  evaluate  the  model 
by  evaluating  how  well  the  fitting  is  and  compute  the  best  fitted  w*,  a, 
and  X  in  the  least-squares  sense  [31]. 

The  procedure  is  described  as  follows:  We  partition  the  data  of 
one  drop  into  two  parts  and  consider  them  separately.  First  we  consider 
only  the  measurements  associated  with  90°  phase  angle.  For  £  =  90°,  we 
have  (from  Eq.  (62)) 

w*2  -  u(90° )2  -  au(90°)3/2  =  0  , 
or  in  another  form, 

a[w(90® )  ] 3/2  =  u*2  -  [u(90® )  ]  2  .  (71) 

By  assuming  an  w*  or,  equivalently,  by  assuming  a  o  (see  Eq.  (16))  and 

plugging  the  measured  w(90°)  into  the  above  equation,  we  obtain  several 

equations  for  one  variable--!*.  We  can  find  an  "optimum"  a  which  fits 

these  equations  the  best  in  the  least-squares  sense. 

The  next  step  is  to  consider  the  measurements  where  £  t  90°.  From 
Eq.  (62),  we  have 

Jfw  -  T(w*2-w2)  ■  a(l+T)w^//^  • 


(72) 
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By  assuming  the  pair  (o,  a)  obtained  in  the  first  step  as  known  (fixed) 
values  and  several  measured  «  at  different  or  the  same  phase  angles,  we 
obtain  several  equations  for  one  variable--*,  whose  optimum  value  is 
then  determined  by  least-squares  method.  Therefore,  for  each  assumed  o, 
we  find  associated  a  and  *.  An  "optimum"  a  or  w*  is  chosen  based  on  the 
overall  performance  of  the  fitting.  The  estimate  for  a  and  *  are 
obtained  along  with  it.  A  general  linear  model  (GLM)  in  a  computer 
system  for  data  analysis  named  SAS  (stands  for  Statistical  Analysis 
System)  supported  by  Yale  Computer  Center  was  used  for  the  least-squares 
regression. 


5.4.2  Examples 

As  an  example,  this  fitting  process  for  the  data  listed  in  Table  2  with 
assumed  o=50.5  dyne/cm  is  shown  in  Appendix  B.  The  result  of  the 
fitting  for  this  set  of  data  is  shown  in  Table  3  and  Fig.  5.3  (Fig.  5.3 
shows  only  the  fitting  for  *).  The  GLM  printout  (Appendix  B)  contains 
some  information  other  than  the  estimated  parameters  (a  and  *) .  The 
information  includes  several  statistical  quantities  that  were  used 
together  to  assess  the  fitting:  the  square  of  the  correlation  function 
(sometimes  called  the  coefficient  of  determination)  (R2),  the 
coefficient  of  variation  (C.V.),  and  the  significance  probability 
<PR>|T|). 

The  square  of  the  correlation  coefficient,  denoted  R2  in  Table  3 
and  R  in  Fig.  5.3,  measures  how  much  variation  in  the  dependent  variable 
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can  be  accounted  for  by  the  model.  The  value  of  the  right  hand  part  of 
Eq.  (71)  or  (72)  obtained  by  employing  experimental  data  is  called  the 
observed  value.  The  value  calculated  by  using  the  fitted  parameters  (a 
or  y)  and  data  to  the  left  hand  part  of  Eq.  (71)  or  (72)  is  called  the 
fitted  or  predicted  value  of  the  dependent  variables.  The  difference 
between  the  observed  and  fitted  value  is  called  the  residual  or  error. 
The  parameter  R2  is  defined  as  the  sum  of  squares  of  the  predicted 
values  divided  by  the  sum  of  squares  of  the  observed  values.  R2  ranges 
from  0  to  1;  the  larger  the  value  of  R2,  the  better  the  model's  fit  (see 
Appendix  B  and  [30]). 

In  Fig.  5.3  and  Table  3  the  symbol  "P"  represents  significance 
probability.  The  significance  probability  is  the  probability  that  the 
observed  data  would  occur  if  a  certain  Null  Hypothesis  were  true;  here, 
the  Null  Hypothesis  is  that  the  parameter  to  be  estimated  (a  or  y)  is 
zero  [32]  [30].  The  significance  probability  ranges  from  0  to  1;  a 
small  value  of  P  indicates  that  the  Null  Hypothesis  is  unlikedly  to  be 
true  or,  in  other  words,  the  model  is  significant. 

The  symbol  CV  in  Table  3  and  the  symbol  c  in  Fig.  5.3  denote  the 
coefficient  of  variation  (C.V.)  which  is  used  to  describe  the  amount  of 
variation  in  the  dependent  variable  (the  right  hand  side  of  Eq.  (71)  or 
(72)).  It  is  equal  to  the  standard  deviation  of  the  dependent  variable 
(after  adjusting  for  the  mean)  divided  by  the  mean  of  the  dependent 
variable  times  100.  The  values  of  C.V.  shown  in  Fig.  5.3  is  normalized 
so  can  fit  into  the  scale  between  0  and  1.  Some  general  observations  of 
the  fitting  are  discussed  below. 


It  can  be  seen  in  Table  3  that  for  a  around  its  "optimum"  value, 
the  fitting  for  o  is  generally  much  better  than  the  fitting  for  y,  which 
should  be  expected  because  the  measurements  for  w(90B)  are  generally 
much  better  than  the  measurements  for  other  phase  angles  in  terms  of 
repeatability.  It  also  shows  that  for  different  assumed  values  of  o  the 
higher  the  a  is,  the  better  the  fitting  for  a  is  ;  on  the  other  hand, 
the  higher  the  a  is,  the  poorer  the  fitting  for  y  is,  and  the  fitting 
for  y  degenerates  first  gradually  then  dramatically  (that  is,  the  R2 
decreases,  while  the  C.V.  and  P  increases  sharply)  as  o  reaches  certain 
value  (Fig.  5.3  shows  an  example). 

Considering  the  fact  that  the  measurements  of  w(90°)  are  the  best 
measurements  we  have,  and  observing  the  fitting  for  both  a  and  y,  we 
choose  the  a  beyond  which  the  fitting  for  y  degenerates  very  quickly  as 
the  optimum  value.  For  most  experiments,  at  the  value  of  "optimum"  a, 
the  fitting  for  y  has  a  R2  about  0.9.  The  value  of  C.V.  and  P  for 
different  test  liquids  could  be  very  different,  because  they  are 
dependent  on  the  mean  of  the  dependent  variable  which  is  a  function  of 
interfacial  tension  and  other  parameters.  However,  the  variation  of 
C.V.  and  P  with  assumed  a  for  individual  experiment  is  clear.  The 
significance  probability  of  the  fitting  for  y  associated  with  the 
optimum  o  is  usually  no  greater  than  0.01,  which  is  the  significance 
level  usually  referred  to  as  "highly  significant"  in  statistics. 


The  optimum  o  we  chose  for  the  above  example  (data  listed  in  Table 
2)  is  50.5,  and  the  associated  a  and  Z  are  1.9  and  13.6  (see  Table  3, 
Fig.  5.3,  and  Table  5).  Employing  the  above  a,  Z,  and  o  obtained  from 
least-squares  fitting,  we  can  reconstruct  the  phase-frequency  relation 
by  the  aid  of  Eq.  (62).  The  result  is  shown  as  the  solid  line  in  Fig. 
5.1.  The  same  procedure  was  applied  to  another  set  of  data  which  is 
shown  in  Fig.  5.2.  The  results  of  the  fitting  for  this  set  of  data  are 
summarized  in  Table  4  and  Fig.  5.4.  The  solid  line  shown  in  Fig.  5.2 
was  reconstructed  based  on  these  "best  fitted"  parameters. 

From  the  results  of  the  fitting  for  the  two  examples  mentioned 
above,  we  may  assess  the  validity  of  the  model  of  Eq.  (62),  especially 
for  describing  the  oscillation  of  the  drop  which  is  formed  in  gel. 

Judged  from  the  values  of  R2Of),  P(JT),  and  C.V.(y)  listed  in  Tables  3 
and  4,  it  appeared  that  the  model  fitted  the  data  well,  even  if  the  test 
drop  was  formed  in  gel  initially.  We  found  that  this  was  generally  true 
for  our  experiments  for  hexane  and  pentane.  The  summary  and  discussion 
of  our  results  for  hexane  and  pentane  are  presented  in  the  next  section. 


TABLE  3 


Fitting  for  one  water-hexane  experiment  (droplet  radius  0.825  mm) 


48.8  48.5  48.8  48.5  58.8  58.5  51.8  51.5  52.8 

INTERFACIAL  TENSION  (DYNE/00 

R  Square  of  the  correction  coefficient 
C  The  coefficient  of  variance  (normalized) 

P  The  significance  probability 

Figure  5.3:  Some  statistics  for  the  least-squares  fitting  for  z-. 

Example  for  host,  water;  drop,  hexane;  drop  radius,  0.0825 
cm  (referred  to  figure  5.1  and  Table  3). 
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TABLE  4 

Fitting  for  one  water-hexane  experiment  (droplet  radius  0.766  mm) 


INTERFACIAL  TENSION  CDYNE/DO 

R  Square  of  the  correction  coefficient 
C  The  coefficient  of  variance  (normalized) 

P  The  significance  probability 

Figure  5.4:  Some  statistics  for  the  least-squares  fitting  for  y : 

Example  for  host,  water;  drop,  hexane;  drop  radius,  0.0766 
cm  (referred  to  figure  5.2  and  Table  4). 


5.5  Measurements  for  Non-superheated  Liquids 
5.5.1  Interfacial  Tension 

A  summary  of  our  interfacial  tension  measurements  for  water-hexane  and 
water-pentane  interfaces  is  presented  in  Table  5.  Parameters  with 

subscript  "th"  are  theoretical  values;  s^  is  calculated  using  Eq.  (68) 

(assuming  the  interfacial  tension  is  known  and  the  value  is  taken  from 

Table  1),  ot^  "is  calculated  by  using  Eq.  (8),  and  y^  is  calcu-  ed  by 

using  Eq.  (9).  Subscript  "f"  denotes  parameters  obtained  usinc  he 

least-squares  method.  sg  is  the  experimental  damping  constant 

calculated  from  data  and  Eq.  (67).  ox  is  the  approximate  interfacial 
tension  calculated  via  Eq.  (70).  o2  is  the  corrected  interfacial 

tension  calculated  using  Eq.  (69)  with  and  y^. 

For  experiments  which  had  the  drops  formed  in  water  (even  the  water 
had  some  gel  dissolved  in  it),  it  appeared  that  the  experimental  damping 

constant  sg  agreed  with  the  theoretical  value  s  ^  well,  also  do  the 

fitted  a,  and  y,  with  a_.  and  y,.  .  For  these  measurements  of  drops 
i  f  th  th 

which  were  first  formed  in  gel,  the  sg  were  much  greater  than  s^;  the 

fitted  ar  is  several  times  greater  than  a.,  while  y,.  is  close  to  y., . 

r  th  r  th 

This  effect  is  expected  because  the  potential  (bulk)  flow  (both  inside 

and  outside  the  drop)  and  the  bulk  viscosity  are  not  changed  much  due  to 

the  presence  of  gel;  on  the  other  hand,  the  flow  in  the  boundary  layer 

is  affected  very  significantly  due  to  some  residue  of  gel  on  the  drop's 

surface  enhancing  the  damping  of  the  drop’s  motion. 
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For  interfacial  tension,  the  corresponded  closely  to  the 

corrected  interfacial  tension  o2.  The  data  in  Table  5  also  demonstrate 
that,  within  the  uncertainty  of  our  measurement,  the  measured 
interfacial  tension  of  the  drops  which  were  first  formed  in  gel  was  not 
different  from  the  measured  value  of  the  drops  which  were  injected 
directly  into  host  water.  Considering  the  similarity  in  chemical 
structure  between  the  superheated  and  non-superheated  sample  liquids 
used  in  our  experiments,  though  without  direct  justification,  we  presume 
that  the  interfacial  tension  between  water  and  the  superheated  liquids 
was  not  significantly  changed  by  the  gel  involved  in  the  transport 
process . 


The  average  measured  interfacial  tension  (o^)  was  50.7  dyne/cm  for 

water-hexane  interface  and  48.6  dyne/cm  for  water-pentane  interface, 
which  were  about  0.8%  lower  than  the  literature  values.  The  total 
uncertainty  of  our  interfacial  tension  measurement  is  estimated  at 
±0.75%.  This  total  uncertainty  is  computed  as  the  square  root  of  the 
sum  of  the  square  of  each  following  error: 

1.  ±0.1%  uncertainty  in  the  frequency  measurement,  which  corresponds 
to  ±0.2%  uncertainty  in  the  measured  interfacial  tension. 

2.  ±0.2%  uncertainty  in  the  droplet  size  measurement,  which 
corresponds  to  ±0.6%  uncertainty  in  the  measured  interfacial 
tension. 

3.  ±0.1%  uncertainty  in  the  density  measurement,  which  corresponds  to 
±0.1%  uncertainty  in  the  measured  interfacial  tension. 
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4.  ±0.4%  uncertainty  in  the  determination  of  the  optimum  a ^  based  on 
the  least-squares  method. 


5.5.2  Damping  and  Viscosity 

In  our  experiments  of  the  quadrupole  oscillation  of  a  drop  in  water,  the 
damping  of  the  motion  is  due  to  the  viscous  dissipation.  Theoretically, 
we  should  be  able  to  deduce  information  on  the  viscosity  of  either  the 
drop  or  the  host  from  the  measured  damping  constant.  Table  6 
illustrates  the  variation  of  the  theoretical  damping  constant 
corresponding  to  several  arbitrarily  assumed  values  of  viscosity.  For 
very  significant  (up  to  30%)  change  in  inner  or  outer  viscosity,  the 

corresponding  variation  of  the  s^  appears  comparable  with  the  typical 

fluctuation  of  the  ratio  s  /s..  (about  15%)  (see  Table  5).  Therefore, 

e  tn 

with  the  precision  of  our  experiment,  it  is  not  practical  to  get  any 
significant  information  about  the  viscosity  of  either  the  inner  or  the 
outer  fluid  from  the  measurements.  The  viscosity  of  0.5%  glycerol  in 
water  solution  is  1%  higher  than  the  viscosity  of  pure  water  [26] .  For 
experiments  marked  with  **  in  Table  5,  the  host  water  had  no  more  than 
0.1%  gel  dissolved  in  it.  From  the  above  example,  the  change  of 
viscosity  of  host  water  due  to  this  amount  of  gel  is  too  small  to  show 
up  in  the  measured  damping  constant,  as  it  is  the  case  we  observed  in 
Table  5. 

For  experiments  denoted  ***  in  Table  5,  if  we  assume  that  near  the 


interface  the  viscosity  was  changed  locally,  then  we  can  estimate  (by 


Eq.  (8))  what  would  be  the  average  viscosity  associated  with  the 

interface  boundary  to  give  th  -  value  a^.  An  example  is  shown  in  Table 

7,  the  value  of  0.03  poise  is  equivalent  to  the  viscosity  of  a  35% 
glycerol  in  water  solution  [26].  The  above  calculation  suggests  the 
existance  of  gel  residue  on  the  interface  for  the  experiments  denoted 
***  in  Table  5,  as  we  discussed  in  section  5.4. 
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TABLE  5 

Measurements  for  non-superheated  liquids. 


nr 

1 

1 

i 

°f 

“f 

*f 

T 

1 

1 

1 

sth 

I 

i 

1 

|water- 

■hexane  interface 

i 

|0.749 

1 

130.63 

1 

-4- 

38.3 

i 

51.3 

50.4 

T 

1 

50.5 

2.3 

16.4 

T 

1 

1 

32.8 

1.8 

17.0 

T 

1 

i 

*  | 

|0.766 

1 

_i_ 

122.14 

1 

_L 

64.5 

51.4 

1 

_L 

KB 

21.3 

1 

| 

31.5  - 

1.8 

16.2 

1 

t 

kkk  { 

1 0.787 

i 

1 

_i_ 

121.63 

1 

1 

1 

33.8 

51.3 

50.5 

1 

_L 

50.5 

2.0 

15.8 

1 

| 

30.0 

1.8 

15.4 

~r 

i 

•k  1 

1 0.825 

i 

1 

_i_ 

113.69 

r 

1 

_L 

30.1 

51.5 

50.6 

l 

_1_ 

50.7 

1.9 

13.6 

1 

1 

27.6 

1.8 

H 

| 

|0.895 

r 

1 

100.23 

1 

1 

28.1 

51.3 

50.5 

_i_ 

50.7 

2.0 

10.1 

1 

1 

23.9 

1.5 

11.9 

~r 

1 

kk  I 

|0.905 

i 

_i_ 

94.01 

~r 

1 

_L 

57.5 

51.5 

■ 

1 

i 

51.1 

4.7 

16.5 

1 

1 

23.4 

1.5 

11.6 

i 

•kkk  I 

|0.985 

r 

1 

_L 

87.20 

-r 

1 

_L 

22.4 

51.4 

50.6 

1 

4. 

50.5 

1.5 

11.8 

T 

1 

_L 

20.1 

1.4 

9.8 

T 

1 

**  1 

j  1.061 

1 _ 

1 

1 

_L 

77.47 

i 

_ L 

23.4 

51.4 

i 

50.7 

r 

1 

JL 

50.7 

1.8 

8.3 

I 

1 

i 

17.6 

1.3 

8.5 

T 

1 

1 

* 

1 

|water-pentane 

i _ 

interface 

jo. 801 

I 

1 

_1_ 

118.08 

i 

1 

i 

30.5 

1 

|  49.4 

48.5 

T 

1 

_L 

48.5 

1.7 

16.5 

"T 

I 

28.1 

1.6 

16.5 

T 

1 

k 

1 0.905 

1 

_j_ 

95.50 

i 

1 

_L 

44.51 

|  49.6 

48.8 

_i_ 

49.0 

3.4 

14.7 

i" 

1 

i 

22.6 

1.4 

12.9 

1 

l 

kkk 

10.710 

1 

1 

_L 

141.18 

1 

1 

_L 

38.9 

|  49.5 

48.5 

7 

_L 

48.5 

2.0 

20.9 

35.0 

1.8 

20.9 

T 

1 

**  i 

| 0.770 

r 

1 

_L 

118.51 

r 

1 

_L 

74.0 

|  49.3 

48.4 

I 

1 

_i_ 

48.7 

5.3 

H 

T 

1 

1 

kkk  1 

J 

1 0 . 791 

1 

1 

_j_ 

120.66 

7 

1 

_L 

29.3 

|  49.5 

48.5 

1 

_L 

48.5 

1.6 

17.5 

i 

28.8 

1.6 

16.9 

1 

1 

1 

*  1 
i 

|0.944 

1 

1 

92.83 

_L 

20.7 

|  49.5 

48.6 

r 

_L 

48.5 

1.2 

12.5 

i 

1 

_L 

20.9 

1.3 

11.8 

1 

1 

i 

kk  j 

1 

1 1.047 

i 

r 

1 

JL 

79.51 

_L 

17.7 

|  49.5 

J _ 

48.7 

_L 

48.7 

1.2 

9.4 

1 

1 

_L 

17.4 

1.2 

9.6 

~T 

_L 

1 

*  ! 

Notes  1)  The  ramark  *  denotes  that  the  drop  was  injected  directly  into  water  from  a 
dispenser  and  there  was  no  gel  dissolved  in  the  water. 

2)  The  remark  **  denotes  that  the  drop  was  injected  directly  into  water  from  a 
dispenser  and  there  was  gel  dissolved  in  the  water. 

3)  The  remark  ***  denotes  that  the  drop  was  formed  first  in  gel  then  transported 
into  water  with  the  gel. 

4)  The  a,  and  have  a  dimension  of  (radian/second)  .  The  y,  and  y_  have  a 

£  tn  £  th 

dimension  of  (radian/second). 
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TABLE  6 

Variation  of  damping  with  viscosity  (I). 


R=0.0825  cm 

a=51.0  dyne/cm 

po=0.9978  g/ml 

Pi=0.6576  g/ml 

yQ=0. 009548  poise 

vi 

poise 

ath 

(rad/sec) 2 

*th 

rad/sec 

Sth 

rad/sec 

0.00313 

1.67 

13.97 

27.64 

0.002 

1.43 

15.52 

25.66  . 

0.004 

1.82 

_ 

13.20 

_ 

28.84 

y^=0. 00313  poise 

Vi 

_ th _ 

Pth 

Sth 

0.009548 

1.67 

13.97 

27.64 

0.007 

1.59 

9.54 

24.49 

0.012 

1.78 

18.53 

30.59 

TABLE  7 

Variation  of  damping  with  viscosity  (II). 


R=0.0766  cm 

Po=0.9973  g/ml  o=51.0  dyne/cm 

Pi=0 . 6576  g/ml 

Po 

ath 

*th 

Sth 

poise 

poise 

is 

(rad/sec) 

rad/sec 

rad/ sec 

0.00313 

0.009548 

1.80 

16.21 

31.59 

0.01 

0.01 

2.60 

14.12 

39.43 

0.02 

_ 

0.02 

_ 

3.68 

28.25 

56.93 

5.6  Measurements  for  Superheated  Liquids 
5.6.1  Test  Liquid 

In  this  section  we  report  measurements  for  butane,  isobutane,  and 
propane.  Some  physical  properties  of  these  three  compounds  are  listed 
in  Table  8.  The  experiments  were  all  conducted  at  the  condition  of  room 
temperature  and  atmospheric  pressure,  and  the  compounds  were  all  in 
liquid  state;  therefore,  they  were  superheated  at  a  degree  from  20°C  to 
over  60aC.  The  boiling  point  at  a  pressure  of  760  mmHg  and  the  vapor 
pressure  at  20°C  are  taken  from  literature  [25].  The  density  in  Table  8 
was  based  on  the  density  at  21°C  and  equilibrium  pressure  [25]  and  then 
corrected  by  using  the  isothermal  compressibility  which  is  extrapolated 
from  the  available  data  at  the  equilibrium  state.  The  surface  tension 
against  its  vapor  listed  in  Table  8  is  taken  from  a  line  which 
extrapolates  the  data  at  the  equilibrium  state  and  connects  the  point  of 
zero  surface  tension  at  the  critical  temperature  (which  should  be  good 
far  from  the  critical  temperature)  (for  example,  Fig.  5.5).  The 
viscosity  listed  in  Table  8  was  extrapolated  from  the  data  at  the 
equilibrium  state.  For  these  values  in  Table  8,  the  uncertainty  is 
estimated  to  be  ±0.2%  for  the  density,  ±1%  for  the  surface  tension  and 
±3%  for  the  viscosity. 
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5.6.2  Measurements  for  Butane  and  Isobutane 

Test  drops  of  butane  and  isobutane  were  formed  in  gel  under  pressure  and 
then  transported  into  water.  The  index  of  refraction  of  butane  (1.3326) 
is  so  close  to  that  of  water  (1.3333)  [26]  that  the  butane  drops  could 
barely  be  seen  by  eyes  while  levitated  in  water,  and  its  shadow  looked 
like  a  ring  with  very  thin  band.  Therefore , very  stable  levitation  and 
fine  adjustment  of  the  stopper  (see  Fig.  3.1)  were  required  to  take 
measurements;  the  signals  from  photomultiplier  for  butane  drops  were 
noisier  than  the  signals  for  drops  of  other  liquids. 

The  results  of  quadrupole  oscillation  experiments  for  butane  and 
isobutane  are  listed  in  Table  9.  As  mentioned  in  the  previous  section, 
for  the  interfacial  tension  deduced  from  data,  ox,  denotes  a  first 

approximation  calculated  using  Eq.  (70),  and  of  is  obtained  using 

least-squares  fittings.  As  shown  in  Table  9,  is  lower  than  Oj  by  an 

average  of  1.6%.  One  can  also  see  that  s  is  two  to  three  times  s..  ,  y,. 

e  th  t 

is  close  to  y^  (within  a  factor  of  two),  and  is  much  greater  than 

atk  (by  a  factor  about  four);  this  behavior  similarly  exhibited  in  our 
data  for  hexane  and  pentane  drops  which  were  formed  in  gel. 

From  our  measurements,  the  average  interfacial  tension  for 
water-butane  interface  is  53.1  dyne/cm,  and  that  for  water-isobutane 
interface  is  54.3  dyne/cm.  Similar  to  the  discussion  in  Section  5.5,  we 
estimate  that  the  uncertainty  for  our  measurements  of  the 
water- isobutane  mterfacia’  tension  is  ±0.8%  with  ±0.2%  uncertainty  in 
the  estimated  density;  the  uncertainty  of  our  measurement  for  the 
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5.6.3  Measurements  for  Propane 

Propane  has  the  highest  degree  of  superheating  among  our  test  liquids. 
Some  special  phenomena  were  observed  in  the  measurements  for  propane: 
the  frequency  associated  with  a  given  phase  angle  increased  with  time, 
while  the  size  of  the  drop  decreased  with  time,  as  can  be  seen  in  Figs. 

5.6,  5.7,  5.9,  5.10,  5.12,  and  5.13.  The  propane  drops  dissolved 
significantly  after  being  levitated  in  water  for  a  couple  of  hours. 

This  fact  makes  it  impractical  to  apply  the  least-squares  method,  as 
described  in  section  5.4,  to  the  data  of  propane  with  water,  because  the 
size  and  frequencies  were  drifting  with  time  and  the  measurements  were 
not  taken  at  the  same  instant. 

Our  approach  for  analyzing  the  data  is  first  to  get  averaged 
variation  of  each  measurement  (e.g.,  the  frequency  associated  with 
certain  phase  or  the  size  of  the  drop)  with  respect  to  time  by  fitting  a 
line  to  the  data  plotted  as  a  function  of  time,  as  shown  in  Figs.  5.6, 

5.7,  5.9,  5.10,  5.12,  and  5.13.  Then  we  can  arbitrarily  choose  an 
instant,  read  the  size  and  frequencies  from  the  fitted  lines,  and 
calculate  the  ox  associated  with  this  particular  instant  (using  Eq. 
(70)).  The  parameters  a  and  JT  can  be  estimated  by  employing  Eqs.  (71) 
and  (72)  with  the  o1(  R,  w(9O0),  and  w(£)  at  a  certain. instant. 

It  was  found  that  the  interfacial  tension  between  water  and  propane 
obtained  in  this  way  was  also  time  depends- :  ;  it  decreased  with  time,  as 

shown  in  Figs.  5.8,  5.11,  and  5.14.  This  decrease  of  interfacial 
tension  with  time  is  presumably  due  to  thr  dissolution  of  propane  in 
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water;  the  concentration  of  propane  in  water  near  the  interface 
increased  with  time,  reducing  the  differential  attraction  force  between 
molecules  across  the  interface. 

The  assumption  that  the  concentration  of  propane  in  the  water  near 
the  interface  increased  with  time  is  supported  by  the  following 
argument:  If  the  diffusivity,  k,  of  propane  into  water  (molecules/unit 
time/unit  area)  is  constant,  then  the  rate  of  change  of  the  drop's 
radius  should  be  inversely  proportional  to  the  cube  of  the  radius,  that 

is,  dR/dt  ~  -k/R^.  But,  as  shown  in  Figs  5.7,  5.10,  and  5.13,  the  slope 
of  the  data  line  is  very  different  from  cubic;  the  rate  of  decrease  of 

the  radius  was  much  less  than  a  cubic  relation  would  predict. 

Therefore,  we  presume  that  the  r.  was  not  a  constant  but  decreased  with 

time,  as  would  occur  if  the  concentration  of  propane  in  water  near  the 

interface  increased. 

Assuming  that  the  true  interfacial  tention  is  1.6%  lower  than  o1( 
as  observed  in  previous  data,  then  the  interfacial  tension  between  water 
and  propane  was  about  58  dyne/cm  when  a  propane  drop  has  contacted  with 
water  for  about  10  minutes  and  continuously  decreased  to  about  54.5 
dyne/cm  when  the  propane  drop  had  remained  in  water  for  two  hours.  We 
estimate  the  uncertainty  of  this  measurement  to  be  ±1%  which  is  of  the 
same  order  as  that  of  measurements  for  butane. 
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Fig.  5.8:  Variation  of  the  interfaoial  tension  between  water  and 
propane  with  for  experiment  #0409. 
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Figure  5.12:  Variation  of  frequency  with  time  for  three  selected 
phase  angles,  experiment  #0408. 


5.6.4  Comparison  with  Fowkes  Theory 

The  results  of  our  measurements  for  interfacial  tension  are  summarized 
in  Fig  5.15.  In  this  figure,  the  square  symbols  represent  data  from  the 
literature;  the  circles  represent  data  from  this  work;  the  curve 
represents  the  calculated  value  given  by  Fowkes1  equation  (referred  to 
Eq.  (1)).  Those  literature  data  of  several  hydrocarbons  with  water 
denoted  by  squares  in  Fig.  5.15  (except  that  of  pentane)  were  used  by 
Fowkes  to  calculate  the  part  of  the  surface  tension  of  water  due  to 
dispersion  force,  which  he  obtained  was  21.8±0.7  dyne/cm  (referred  to 
section  1.6  and  [13]).  The  surface  tension  of  the  superheated  liquids 
are  taken  from  Table  8.  The  interfacial  tension  between  water  and 
propane  shown  in  this  figure  is  58  dyne/cm  which  corresponds  to  the 
measurements  taken  at  an  "early  stage"  in  the  period  of  measurement 
(refer  to  section  5.6.3).  Fig.  5.15  demonstrates  that  for  the  three 
superheated  liquids,  the  agreement  between  our  data  and  Fowkes1  theory 
is  fair,  similar  to  that  between  the  data  of  pentane  with  water  and 
Fowkes'  prediction;  the  trend  of  increasing  interfacial  tension  with 
decreasing  surface  tension  in  this  region  of  lower  surface  tension  is 
clearly  shown  in  both  data  and  Fowkes'  theory. 
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Figure  5.t5:  Comparison  between  dota  and  Fowkes  theory, 
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5.7  Measurements  for  Compressadensity 

Applying  the  established  acoustic  levitation  technique,  we  measured  the 
compressadensity,  which  is  a  bulk  property,  of  the  three  superheated 
liquids.  The  procedure  is  outlined  in  section  4.5,  and  the  calculation 
is  based  on  Eqs.  (60)  and  (61). 


5.7.1  Control  Experiment 


As  a  test  of  the  accuracy  of  the  levitation  technique  for  measuring 
adiabatic  compressibility  (g)  (assuming  the  density  is  known),  control 
experiments  were  performed  in  which  one  of  the  two  liquids,  n-heptane 
and  n-pentane,  was  considered  an  unknown,  and  the  other  liquid  a 
reference.  The  sound  speed  (C)  in  the  "unknown"  liquid  was  measured  by 
a  direct  method  (see  section  3.5).  Making  use  of  the  relation  C2  = 
l/pf5,  we  can  compare  the  results  obtained  from  levitation  technique  with 
that  obtained  from  the  direct  method.  The  results  are  summarized  in 
Table  10.  The  densities  of  n-heptane  and  n-pentane  in  Table  10  were 
measured  using  the  Mettler/Parr  density  meter.  The  density  of  water  is 
taken  from  the  literature  (e.g.  [26]).  The  sound  speeds  in  the  three 
liquids  were  measured  using  Nusonics  sonic  solution  monitor  (see  section 
3.5). 


Our  data  show  that  the  precision  of  our  measurements  for 
compressibility  is  about  1%,  and  for  sound  speed  is  about  0.5%.  The 
overall  uncertainty  of  this  measurement  is  composed  of  the  following 
errors  [5] : 


1.  ±0.2%  uncertainty  in  the  voltage  measurement,  which  is 
attributable  to  the  precision  of  the  digital  ac  voltmeter  and  the 
deviation  of  the  actual  experimental  condition  from  the  assumption 
that  the  pressure  of  the  incident  waves  for  both  the  unknown  and 
reference  drops  is  proportional  to  the  input  voltage  to  the 

t "  .nsducer . 

2.  ±0. 002cm  uncertainty  in  droplet  position  monitoring,  which 
corresponds  to  ±0.2%  uncertainty  in  voltage  measurement. 

3.  ±0.1%  uncertainty  in  density  measurement. 

4.  ±0.1%  uncertainty  in  sound  speed  measurement  for  host  and 
reference  liquid. 

The  compressibility  and  sound  speed  measured  using  the  levitation 
technique  agreed  well  with  that  obtained  using  traditional  technique, 
with  the  difference  within  the  uncertainty  of  our  levitation 
measurements . 


5.7.2  Compressadensity  of  Superheated  Liquids 

The  same  procedure  was  applied  to  butane,  isobutane,  and  propane  liquid, 
which  were  superheated,  to  measure  their  adiabatic  compressibility  and 
sound  speed.  Two  references,  n-heptane  and  n-pentane,  were  used  in  all 
measurements.  The  densities  of  the  superheated  liquids  are  taken  from 
Table  8.  The  results  are  listed  in  Table  10.  We  estimate  the  total 
uncertainty  to  be  about  1%  for  sound  speed  and  2%  for  compressibility. 
The  results  shown  in  Table  10  demonstrate  that  the  results  for  the  two 
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different  references  are  consistent  with  each  other,  the  difference 
being  within  the  uncertainty  of  this  measurement.  Since  the 
compressibility  and  sound  speed  are  bulk  properties  of  liquids,  we 
judged  that  the  data  should  not  be  affected  by  the  presence  of  little 
amound  of  residuce  of  gel  on  the  interface.  There  are  no  other  data 
available  for  the  results  of  the  superheated  liquids  listed  in  Table  10 
to  compare  with.  However,  the  trend  that  the  liquid  of  higher  degree 
superheating  behaves  closer  to  being  in  gaseous  state  (greater 
compressibility  and  lower  sound  speed)  is  clear. 


Chapter  6 


SUMMARY  AND  CONCLUSIONS 

In  this  work,  we  have  studied  a  technique  for  measuring  the  interfacial 
tension  between  two  liquids  and  applied  it  to  the  interfaces  between 
water  and  each  of  three  superheated  liquids.  The  method  of  measurement 
is  based  on  the  resonance  properties  of  a  drop  which  is  levitated  in 
another  liquid  (the  host)  and  also  driven  into  quadrupole  shape 
oscillation  by  an  acoustic  lorce.  The  shape  oscillation  technique 
described  here  has  the  special  advantage  that  the  host  liquid  provides  a 
very  smooth  and  clean  container  for  the  levitated  drop,  thereby  greatly 
improving  its  time  of  survival  in  the  superheated  state. 

We  have  developed  a  simplified  model  to  describe  the  motion  of  a 
drop  oscillating  in  a  quadrupole  mode  in  a  host  liquid.  In  this  model, 
the  flow  field  is  partitioned  into  several  parts:  the  potential  (bulk) 
flow  in  the  outer  (host)  and  inner  (drop)  fluid  as  well  as  the  boundary 
layer  in  the  outer  and  inner  fluid  contiguous  to  the  interface;  they  are 
superimposed  in  an  ad  hoc  manner  to  describe  the  whole  flow.  The 
phase-frequency  relation  deduced  from  this  simplified  model  agrees  well 
with  that  given  by  previous  theories  (e.g.  Marston's)  which  were  derived 
in  an  exact  approach.  This  simplified  model  reveals  the  physical  nature 
of  the  problem  more  explicitly.  Though  the  derivation  of  this 
simplified  model  was  specified  to  the  quadrupole  mode,  the  procedure 
should  be  applicable  to  oscillations  of  higher  modes. 
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The  apparatus  of  this  experiment  is  similar  to  that  of  some 
previous  work,  with  some  modifications.  The  accuracy  of  the  size 
measurement  has  been  improved  and  it  contributes  significantly  to  the 
precision  of  our  measurements  for  interfacial  tension.  Our  results  for 
the  interfaces  between  water  and  each  of  two  common  liquids  were  lower 
than  the  literature  values  by  about  0.8%.  The  uncertainty  of  the 
measurements  was  estimated  to  be  ±0.75%.  Though  the  damping  of  the 
drop's  oscillation  is  due  to  viscous  dissipation,  it  appeared  that, 
within  the  precision  of  our  experiment,  it  is  impractical  to  get  any 
significant  information  about  the  viscosity  of  either  the  drop  or  the 
host  liquid  from  the  measurements. 

We  developed  a  process  for  introducing  a  liquid  drop  into  water  by 
the  aid  of  a  gel  which  is  immiscible  with  the  drop  but  soluble  in  water. 
It  is  a  convenient  way  to  handle  superheated  liquids.  However,  test 
drops  that  underwent  this  process  exhibited  much  greater  damping  in 
their  shape  oscillation  and  this  fact  raised  questions  as  to  whether  we 
could  infer  the  interfacial  tension  from  measurements  by  using  the 
approximation  method  used  in  a  previous  work.  We  therefore  developed  a 
method  based  on  the  least-squares  principle  for  estimating  interfacial 
tension  and  damping  parameters.  The  results  of  this  least-squares 
method  for  two  common  liquids  demonstrate  that  the  interfacial  tension 
and  the  viscous  dissipation  in  the  bulk  flow  were  not  affected  by 
detectable  amounts  due  to  the  presence  of  the  gel,  while  the  viscous 
dissipation  in  the  boundary  layer  increased  very  significantly. 


Measurements  of  interfacial  tension  were  reported  for  the 
interfaces  between  water  and  butane,  isobutane,  as  well  as  propane. 
These  three  compounds  were  in  the  superheated  liquid  state  under  our 
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experimental  condition  (room  temperature  and  atmospheric  pressure).  The 
uncertainty  of  the  measurements  was  estimated  to  be  ±1%.  The  size  of 
propane  drops  decreased  gradually  over  a  two  hour  period  in  water;  the 
interfacial  tension  decreased  with  time  due  to  the  increase  of  propane 
concentration  in  water  near  the  interface.  Our  results  show  fair 
agreement  with  the  prediction  given  by  Fowkes1  theory. 

We  have  reported  the  sound  speed  and  compressibility  of  those  three 
superheated  liquids  which  were  measured  using  an  established  acoustic 
levitation  technique.  The  uncertainty  is  estimated  at  ±1%  for  the  sound 
speed  measurement  and  at  ±2%  for  the  compressibility  measurement.  Our 
results  show  that  the  sound  speed  decreases  and  the  compressibility 
increases  as  the  degree  of  superheat  increases,  as  one  should  expect. 

This  work,  as  one  among  others,  demonstrates  that  the  acoustic 
levitation  technique  is  a  useful  tool  for  studying  either  surface 
properties  (e.g.  interfacial  tension)  or  bulk  properties  (e.g.  sound 
speed  and  compressibility)  of  fluids  with  a  special  advantage  in  dealing 
with  metastable  liquids.  This  work  reports  data  v/hich  were  previously 
unavailable  for  the  interfacial  tension  and  compressibility  of  three 
superheated  liquids  and  may  have  applications  in  various  fields 
(referred  to  page  1,  5,  and  11).  Though  we  judged  that  the  measured 
interfacial  tensions  were  not  affected  by  the  presence  of  gel,  any 
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uncertainty  on  this  point  can  only  be  removed  by  conducting  experiments 
without  using  gel  or  other  material  in  the  preparatory  steps  prior  to 
the  actual  interfacial  tension  measurements. 

Some  future  work  relating  to  this  technique  is  desirable.  The 
approach  of  the  simplified  analytical  model  may  be  extented  to  some 
special  situations  in  which  the  exact  approach  may  have  difficulty  due 
to  the  geometry,  such  as  when  the  equilibrium  shape  of  the  drop  is 
non-spherical .  It  is  also  possible  to  superimpose  terms  such  as 
interfacial  viscosity  and  interfacial  elasticity  into  this  simplified 
model  for  describing  some  situations  such  as  when  there  is  a 
surface-active  film  or  membrane  between  the  two  fluids;  with  this 
modification  in  the  analytical  model  and  by  performing  measurements  as 
described  in  this  manuscript,  we  should  have  a  sensitive  method  of 
characterizing  the  surface  layer. 


APPENDIX  A 


Methods  for  Data  Reduction 

The  model  that  we  used  to  describe  the  motion  of  an  oscillating  drop  in 
a  liquid  host  is  expressed  in  Eq.  (62), 

tan£  =  (aw^^  +  Jw)/(u)*2  "  -  u2)  ,  (62) 

where  £  and  «  are  measured  in  our  experiments  (see  section  5.2.1).  By 

measuring  three  frequencies  (w)  associated  with  three  different  phase 

angles  (£),  and  by  substituting,  them  into  the  above  equation,  we  get  a 

system  of  three  linear  equations  and,  theoretically,  we  can  solve  the 

three  unknowns:  a,  2f,  and  w*2  (or  a)  by  direct  methods  such  as  Cramer's 

rule  or  Gauss  elimination.  However,  it  appears  that  the  system  of 

equations  mentioned  above  is  “ill-conditioned",  and  the  direct  methods 

are  of  very  little  use  for  inferring  interfacial  tension  from  the 

phase-frequency  measurements  with  the  precision  we  can  achieve  in  our 

measurements.  An  example  is  shown  below  in  two  steps: 

(i).  Making  use  of  Eq.  (16),  we  may  rewrite  Eq.(62)  in  the 
following  form: 

[ (T«24)/ (R3T) ] a  -  («)r  -  [(l+T)w3/2]a  =  Tw2  .  (73) 

Assuming  R=0.0825  cm,  ct=1.88,  y=13.56,  a=50.7  (dyne/cm)  (see  Table  5), 

Pq=0.9978,  p^=0.6576  (see  Table  1),  and  employing  Eq.  (62),  we  may  get  a 

set  of  "artificial  data",  such  as: 

f(Hz)  s(deg.) 


111.729 

114.081 

117.609 


68.257 

95.180 

128.363 
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Substituting  this  set  of  artificial  data  into  Eq.  (73),  we  get  a  system 
of  linear  equations  [A]  [X]  =  [B] ,  where 


and 


[X]  =  [o  y  a]T 

r  27005.89 

-702.02 

[A]  =  j -118809 . 38 

-716.79 

L  -13606.97 

-738.95 

r  1235708.00-, 

[B]  =  | -5667649 . 00 | 

<-  -689871. 87-1 

-65238.53-, 
192502.12 | 
5290.33-1 


Employing  Cramer's  rule,  we  solve  the  unknown  matrix, 

iT 


[X]  =  [50.71  13.50  1.91] 

which  is  very  close  to  the  original  assumed  values. 


(74) 


(ii)  If  we  make  a  small  change  in  the  artificial  data,  say,  one  of 
the  frequency  from  117.609  to  117.800,  then  it  produces  "small"  changes 
in  the  coefficients: 

r  27005.89  -702.01  -65238.53-, 

[A]  =  j - 118809 . 37  -716.79  192502. 12  j 

L  -13606.97  -740.15  5303. 22-1 

r  1235708. 00-, 

[B]  =  j -5667649 . 00  j 

L  -692114.56-1 

But  the  corresponding  changes  of  the  solution  are  very  "large". 

[X]  =  [60.67  -126.41  7.54]T. 

The  performance  of  direct  methods  can  be  predicted  by  a  theory  in 
numerical  analysis  (e.g.  [33])  which  provides  an  estimate  of  the 

relative  error: 

1  |AXc  -  B  ||  ||Xc  -  Xt  ||  ||AXc  -  B  || 

- - - <  -  <  || A  ||  || A~ 1  ||  - 

II ^  II  II A  ||  || B  ||  ||Xt||  || B  || 


where  1*11  denotes  a  matrix  norm,  subscript  c  denotes  the  computed  value, 
subscript  t  denotes  the  true  value,  and  all  variables  are  matrix.  The 

term  ||A||«||A  *||  is  usually  called  the  condition  number  and  its  size  is  a 
measure  of  how  good  a  solution  we  can  expect  direct  methods  to  generate. 

The  condition  number  satisfies  the  inequality  ||A||«||A  1 1|  i  1.  Following 
a  specific  definition  of  the  matrix  norm[33] ,  we  can  calculate  the 

condition  number  for  the  above  example  (referred  to  Eq.  (74)),  such  as: 

»A|t  •  HA*1  II i  =  1.6  x  10*  »  1 
||A J  .  j|A‘Y  =  3.6  x  10*  >>  1 

The  above  calculation  indicates  that,  with  the  precision  we  can  achieve 
in  our  measurements,  if  we  use  direct  methods  to  solve  [X],  the  relative 
error  could  be  unacceptably  large. 

The  problem  can  also  be  demonstrated  by  the  following  observations. 
First,  the  typical  magnitudes  of  the  terms  in  Eq.  (62)  are  very 
different,  as  shown  in  the  following  vector  diagram  which  is  drawn 

approximately  to  scale  for  the  above  example  with  f  =  114.081  H  and  £  = 

z 

95.180  degree. 


a(w3/2)/2 

,  * 

W2  u*2 

where  w*2  =  546062,  w2  =  513791,  yu  *  9720,  and  a(w3/2)/2  =  51023. 

3/2 

Secondly,  the  dependence  on  w  of  the  two  damping  terms  (Jw  and  aw  )  is 
rather  close,  and  the  variation  in  frequency  is  rather  small  for  the 

range  of  phase  angles  covered  in  our  measurements. 


The  idea  of  our  approach  to  analyze  the  data  is  formed  by  the 
following  reasoning.  To  avoid  the  difficulty  in  direct  methods  and  to 
use  all  the  data  we  have  for  one  experiment  (usually  more  than  three 
phase-frequency  measurements),  we  employ  a  least-squares  principle  to 
fit  the  model  in  Eq.  (62)  to  our  data.  To  give  greater  emphasis 

3/2 

relatively  to  the  terms  Xu  and  ct(w  ),  we  assume  that  u>*  is  known  by 
guessing  a  value  for  the  interfacial  tension.  Then  the  terms  w*2  and  w2 

almost  cancel  out  each  other.  In  the  remaining  two  unknowns,  X  and  a, 

only  a  affects  the  value  of  f(90°)  (see  Eq.  (71));  also,  our 

measurements  for  f(90°)  are  usually  much  better  than  f(£),  where  £*90°, 

in  terms  of  repeatability.  Therefore,  we  first  employ  the  least-squares 

method  and  only  the  data  of  f(90°)  to  get  the  best  fitted  «.  We  then 

use  the  assumed  u»*,  fitted  a,  and  the  data  of  f(£*90°)  to  obtain  the 

best  fitted  y.  The  assumed  «*  which  is  the  closest  to  the  "true"  value 

should  lead  to  the  best  fitting  (referred  to  section  5.4).  Details  of 

the  procedure  are  described  in  section  5.4  and  Appendix  B. 
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APPENDIX  B 

An  Example  for  the  Least-squares  Fitting 

This  appendix  provides  an  example  for  the  process  of  least-squares 
fitting  that  we  employed  to  analyze  our  data.  The  first  part  of  the 
attached  printout  is  the  data  set  used  in  this  example  (see  Table  2). 

The  remaining  part  of  this  printout  is  the  GLM  printout  (see  section 
5.4).  Relevant  terms  are  marked  with  numbers  and  are  explained  as 
follows : 

(1) .  the  data  used  in  this  example  (see  Table  2) 

(2) .  the  assumed  interfacial  tension. 

(3) .  the  coefficient  of  a  in  the  left  hand  part  of  Eq.  (71). 

(4) .  the  right  hand  part  of  Eq.  (71). 

(5) .  the  sum  of  squares  for  the  dependent  variable  in  (4). 

(6) .  the  sum  of  squares  for  the  predicted  dependent  variable  in  (10). 

(7) .  R2,  the  ratio  produced  by  dividing  (6)  by  (5). 

(8) .  the  difference  between  (4)  and  (10). 

(9) .  the  sum  of  squares  of  (8). 

(10) .  the  value  obtained  by  employing  the  estimated  parameter  (15)  into 

Eq.  (71). 

(11) .  (9)  divided  by  the  DF  (degree  of  freedom). 

(12) .  the  square  root  of  (11). 

(13) .  the  mean  value  of  (4). 

(14) .  (12)  divided  by  (13),  then  times  100,  (the  coefficient  of  variation). 

(15) .  the  estimated  parameter  based  on  minimizing  (9). 

(16) .  the  value  for  testing  the  null  hypothesis  that  the  parameter 

equals  zero,  which  equals 


(13)-0.0 

(12)//n 


where  n  is  the  number  of  observations  (n=3  in  this  example)  [32]. 

« 

(17).  This  value  answers  the  questions,  "If  the  parameter  is  really 

equal  to  zero,  what  is  the  probability  of  getting  a  larger  value 
of  (16)?"  Based  on  the  assumption  that  the  dependent  variable  has 
a  normal  distribution  with  mean  zero  and  estimated  standard 
deviation  obtained  from  n  observations  (the  term  (12)),  the  value 
of  (17)  is  then  the  probability  to  get  a  greater  value  of  (16)  or, 
equivalently,  a  greater  value  of  (13)  [30]  [32]. 
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*****5ubstitute  the  data  in  Table  2  into  Eq.  (71), 

assume  the  interfacial  tension  to  be  50.5  dyne/cm,  and 
invoke  SAS  GLM  program  to  fit  alfa:  |2) 

STATISTICAL  ANALYSIS  SYSTEM 

OBS  ALFA  INER 

(3)  (4) 

1  190.795  340.816 

2  190.946  335.434 

3  190.996  333.641 

STATISTICAL  ANALYSIS  SYSTEM 

GENERAL  LINEAR  MODELS  PROCEDURE 
DEPENDENT  VARIABLE:  INER 


SOURCE 

DF 

SUM  OF  SQUARES 

ii 

MEAN  SQUARE 

F  VALUE 

MODEL 

1 

'6* 339956.84381956 

339956.84381956 

22129.25 

ERROR 

2 

(®)  30.72466672 

^^15.36233336 

PR  >  F 

UNCORRECTED 

TOTAL  3 

.  .339987.56848628 

(5) 

0.0001 

R-SQUARE 

tn\ 

(u)  c  v- 

(12)  STDDEV 

INER  MEAN 

(is) 

0.999910'  ' 

1.1643 

3.91948126 

336.63020000 

SOURCE 

DF 

TYPE  I  SS 

F  VALUE  PR  >  F 

ALFA 

1 

339956.84381956 

22129.25  0.0001 

S  T  A  T  I  S 

T  I  C  A  L  AN 

ALYSIS  SYS 

T  E  M 

GENERAL  LINEAR  MODELS  PROCEDURE 
DEPENDENT  VARIABLE:  INER 


SOURCE 

ALFA 


TYPE  IV  SS  F  VALUE 


1  339956.84381956  22129.25 


PR  >  F 

0.0001 


PARAMETER 

ALFA 


ESTIMATE 


1.76326420 


T  FOR  HO: 
PARAMETERS 

(16) 

148.76 


PR  >  | T | 

(17) 

0.0001 


STD  ERROR  OF 
ESTIMATE 

0.01185316 


OBSERVATION 

OBSERVED 

PREDICTED 

.  ,  RESIDUAL 

VALUE 

(10)  VALUE 

(8) 

1 

340.81640000 

336.42146348 

4.39493652 

2 

335.43360000 

*  336.68806903 

-1.25446903 

3 

333.64060000 

336.77676122 

-3.13616122 

STATISTI 

CAL  ANALYS 

IS  SYSTEM 

GENERAL  LINEAR  MODELS  PROCEDURE 
DEPENDENT  VARIABLE:  INER 

SUM  OF  RESIDUALS  '  0.00430626 
SUM  OF  SQUARED  RESIDUALS  30.72466672 
SUM  OF  SQUARED  RESIDUALS  -  ERROR  SS  -0.00000000 
FIRST  ORDER  AUTOCORRELATION  -0.05139501 
DURBIN- WATSON  D  1.15400921 


*****  The  end  of  SAS  GLM  output 
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